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ABSTRACT 



Context. Imaging surveys of dust continuum emission at (sub)millimetre wavelengths provide a powerful tool to study molecular 
clouds and the early stages of star formation. 

Aims. Through submm dust continuum mapping, we attempt to search for genuine infrared dark clouds (IRDCs) and precursors to 
massive stars and stellar clusters in the Galactic plane, and to determine their basic physical properties. 

Methods. We have mapped four selected about 0°5 x ffi5-sized fields containing Spitzer 8-pm dark regions with APEX/LABOCA at 
870 fim. Selected positions in the fields were observed in C"0(2 - 1) to obtain kinematic information. The obtained LABOCA maps 
are used in conjunction with the Spitzer IR images. 

Results. The total number of clumps identified in this survey is 91, out of which 40 (44%) appear dark at 8 and 24 pm. The remaining 
clumps are associated with mid-IR emission. Seven clumps associated with extended-like 4.5 pm emission are candidate extended 
green objects (EGOs). Filamentary dust "ridges" were found towards the Spitzer bubbles NlO/11 in one of our fields. The relative 
number of IR-dark and IR-bright clumps suggest that the duration of the former stage is about 1.6 x 10' yr The mass distribution 
of the total sample of clumps, and that separately constructed for the IR-dark and IR-bright clumps, could be fitted at the high-mass 
end with the power-law function dN/d log M oc M"^, where F ^ 0.7 ... 0.8. The C"0 observation positions appear to be dominated 
by non-thermal motions, and the data also revealed some potential sites of strong CO depletion. In Gl 1.36+0.80, which is the best 
example of a filamentary IRDC in our sample, the clumps appear to be gravitationally bound. The fragmentation of the filament can 
be understood in terms of a "sausage"-type fluid instability, in agreement with the results for other IRDCs. The fragmentation and the 
CO depletion timescales in Gl 1.36 appear to be very similar to each other 

Conclusions. Many of the identified clumps are massive enough to allow high-mass star fomation, and some of them already show 
clear signposts of that. In the NlO/11 bubble environment, the morphology of the detected dust emission conforms to the triggered 
high-mass star formation in the system. The clump mass distributions are similar to those found for diffuse CO clumps, and can be 
explained by the action of supersonic turbulence. The formation of filamentary IRDCs might be caused by converging turbulent flows, 
and the same process may play a role in exciting the fluid perturbations responsible for the fragmentation of the clouds into clumps. 

Key words. Stars; formation - ISM: clouds - Submillimetre: ISM 



1. Introduction 

Understanding the origin of high-mass (M > 8 M0) 
stars, particularly the first steps in the formation process, 
is probably one of the greatest challenges of modern astro- 
physics. Since the discovery of the so-called infrared dark 
clouds, or IRDCs dPerault et al. 19961 [Eganetal. 1998) 1, am- 
ple evidence has been gathered concerning their important 
role in the earliest stages of Galactic high-mass star for- 
mation (e.g.,^athbome et al. 2006, Beuther & Steinacker 2007] 
[Chambers et al. 2009. Battersby et al. 2010, ,Zhang et al. 20 iT] 
and many other works). In particular, studies of IRDCs have the 
potential to help to understand the inital conditions of high-mass 
star and stellar cluster formation, which is necessary in order 
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* This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX) under programme 087.F-9315(A,B). 
APEX is a collaboration between the Max-Planck-Institut fiir 
Radioastronomie, the European Southern Observatory, and the Onsala 
Space Observatory. 



to constrain, or even distinguish, between different theoretical 
view£l 

Thermal dust continuum emission at far-infrared (FIR) 
and (sub)millimetre wavelengths provides a powerful observa- 
tional tool to search and study the densest parts of IRDCs. 
Dust continuum imaging with bolometer cameras can be used 
to distinguish the real IRDCs, i.e., cold dense molecular 
clouds, from the minima in the Galactic mid-infrared (MIR) 
background radiation, which may look like candidate IRDCs 
(IWilcock et al . 2012'). Optically thin dust emission also provides 
a probe of the basic physical properties of dense clouds, such as 
the column density of molecular hydrogen and the mass of the 
cloud. This information is needed to learn the physical condi- 
tions prevailing in the precursor regions of stellar clusters and 
high-mass stars. 

In this paper, we present the results of our submm dust 
continuum observations at 870 fim of four selected regions in 
the Galactic plane, each of which contain IRDCs. Throughout 



' At least in some high-mass star-forming regions, however, other 
factors, such as interactions between the cluster members, stellar feed- 
back, and external forces, may be more important for the source evolu- 
tion than the initial conditions of the parent cloud. 
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the paper, we will use the term "clump" to refer to sources 
whose typical radii, masses, and mean densities are, respec- 
tively, - 0.2 -1 pc, ~ 10^ - 10^ Mq, and 10^ - 10* cm^^ 
(cf. |Bergin & Tafalla 2 007 ). The rest of the present paper is 
organised as follows. Observations and data reduction are de- 
scribed in Sect. 2. Observational results are presented in Sect. 3. 
Analysis and its results are presented in Sect. 4, and in Sect. 5 we 
discuss the obtained results. In Sect. 6, we summarise the results 
and draw our main conclusions. 



2. Observations and data reduction 

2.1. Archival data from the Spitzer Space Telescope 



In this study, we use the Spitzer dWerner et a l. 2004] l IR 
data taken as part of the GLIMPSE (Benjami n et al. 2003 
IChurchwell et al. 20 09^ and MIPSGAL (Care y etal. 2009 
Galactic plane surveys. The former survey employed the 
IRAC instrument operating at 3.6, 4.5, 5.8, and 8.0 yum 
dFazio et al. 2004l l. whereas the latter one used the MIPS 
instrument at 24 and 70 yum (IRieke et al. 20041) . The angular 
resolution of the Spitzer-IRAC instrument is V.'9 at 8 //m, and 
that of the MIPS instrument is 6" at 24 fim. We note that we 
have used the data provided by both the GLIMPSE I and II 
surveys, which covered the nominal Galactic longitude ranges 
of 10° < |/| < 65° and \l\ < 10°, respectively. The data were 
retrieved from the Spitzer science archiv^ 

The 8 yum images of the Galactic plane are particularly useful 
for the search of candidate IRDCs. The GLIMPSE 8 //m band 
contains the UV-excited 7.7 and 8.6 fim PAH (poly cyclic aro- 
matic hydrocarbon) features (e.g.. lDraine 20 03 1. which together 
with emission from warm interstellar dust yield a bright MIR 
background. High columns of cold dust in IRDCs cause them to 
appear as dark absorption features against this background radi- 
ation field. 

2.2. LABOCA dust continuum mapping 

As a starting point of our study we visually inspected the 
Spitzer-GLIMPSE 8-jum images of the Galactic plane, and chose 
four target fields containing filamentary IR-dark features to be 
mapped in the submm dust continuum emission. Because IRDCs 
often exhibit filamentary shapes and are relatively devoid of (vis- 
ible) star formation, the target sources of this study are likely 
to represent rather typical IRDCs. The selected fields, which 
all belong to the first Galactic quadrant (0° < I < 90°), are 
listed in Table[T] These target fields were mapped with the Large 
APEX Bolometer CAmera (LABOCA; [Siringo et al. 2009| l on 
the 12-m Atacama Pathfinder Experiment (APEX) telescope 
at Llano de Chajnantor in the Atacama desert of the Chilean 
Andes ( Glisten etal. 2006fl The LABOCA instrument is a 
multi-channel bolometer array, where 295 semiconducting com- 
posite bolometers are arranged in a series of nine concentric 
hexagons around a central channel. The system operates at a 



^ 'http : //irsa . ipac . caltech . edu/data/SPITZER/docs/ 1 
spitzerdataarchives/ 

^ The observations presented here are not part of the APEX 
Telescope Large Area Survey of the GALaxy (ATLASGAL) 
conducted with the LABOCA array (Schull er et al.'2069] 
http://www.mpifr.de/div/atlasgal/index.html ). However, 
the ATLASGAL survey covers the same regions we have mapped. At 
the time of writing, the ATLASGAL data were not yet released into the 
public domain. We note that ATLASGAL has the average Icr rms of 
~ 50 mJy beam"' and the resolution of ~ 19'.'2. 



central frequency of 345 GHz (A - 870 yum) with a bandwidth of 
about 60 GHz to match the coiTesponding atmospheric window. 
The nominal angular resolution of LABOCA is 19'.'2 ± a.'3 (half 
power beamwidth; HPBW), and its total field of view is 1 1 '4 
(about 0.09 pc and 3.3 pc at 1 kpc, respectively). 

Our LABOCA observations took place on 1 9 May 20 1 1 , dur- 
ing the UTC time ranges of 03:17-06:20 and 07:54-10:50. The 
observing conditions were very good: the atmospheric zenith 
opacity, as determined using skydip measurements, was in the 
range = 0.10 - 0.12, and the amount of precipitable water 
vapour (PWV) was in the range 0.15-0.30 mm. The telescope 
focus and pointing were, respectively, optimised and checked at 
regular intervals on the planets Saturn and Neptune, the Class 
protostellar core IRAS 16293-2422, the ultracompact (UC) 
Hii region GIO.620.38, and the massive young stellar object 
(MYSO) G305. 80-0.24 (B13134). The absolute calibration un- 
certainty is estimated to be about 10%. 

The observations were performed using the on-the-fly (OTF) 
mapping mode, in which the telescope scanned continuously in 
right ascension (RA) along each row. We used a scanning speed 
of 3' s"' and the step size 6'.'5 (~ 1/3 the beam HPBW) be- 
tween RA subscans. The step size < l/3x the beam HPBW 
is recommended to avoid beam broadening. The angular sizes 
of the maps are given in Col. (2) of Table [1] (sizes are in the 
range of ~ 0.23 - 0.30 deg^, with a total angular area of about 1 
deg^). The target fields were mapped three to five times, with to- 
tal on-source integration times in the range 49-79 min [Col. (4) 
of Table 

Data reduction was done with the CRUSH-2 
(Comprehensive Reduction Utility for SHARC-2) (version 
2.11-al) software packag^ (IKovacs 20081) . We used the 
pipeline iterations both with the default reduction parameters, 
and also with specifying the 'extended' option, which better 
preserves the extended structures. In the cases of Gl. 87-0. 14, 
G2.llH-0.00, and G13.22-0.06 (hereafter, G1.87, etc.), the 
'extended' -reduced maps were finally chosen for the analysis 
because fainter extended structures were clearly better recov- 
ered compared to the default reduction method. In the case 
of G11.36, however, we adopted the map reduced with the 
default parameters because the clumpy structure of the filament 
became more clearly visible (Fig. |3}- A slight beam-smoothing 
was applied in the reduction process, i.e., the maps were 
smoothed with a Gaussian kernel of the size 3'.'8 (full width at 
half maximum; FWHM). The instrument beam HPBW used 
by CRUSH-2 was 19'.'5, and therefore the angular resolution 
of the final maps is 19'.'9 (~ 0.1 pc at 1 kpc). The gridding 
was done with a cell size of 4". The resulting Icr rms noise 
levels in the final co-added maps are ~ 40 - 90 mJy beam ' 
[Col. (5) of Table [T|. Assuming that the 870-yum dust opacity 
and the dust temperature are, respectively, 1.38 cm^ g"' and 15 
K (see Sect. 4), the above surface-brightness sensitivity levels 
translate into Icr H2 column-density detection thresholds of 
^(Ha) - 2.0 - 4.4 x 10^' cm"^. These correspond to visual 
extinction values of Ay = A^(H2)/0.94 x 10^' - 2.1 - 4.7 mag 
(IBohlin et al. I978lFl We note that employing the 'extended' 
option in the reduction process leads to maps with a higher noise 
level than the values of ~ 30 - 50 mJy beam"' resulting from 
the standard procedure (because large scale emission is tried to 
be preserved). Therefore, the noise in the map of Gl 1.36, which 



http : //www . submm . caltech . edu/~ share/crush/ index . htmj 
' This Av - ^(112) relationship is based on observations of diffuse 
interstellar medium, and therefore may not be exactly correct for dense 
molecular clouds. 
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Table 1. Target fields mapped with LABOCA. 



Field" 


Map size 


No. of maps 


On-source time 


lo" rms noise*" 




['X'] 




[min] 


[mJy beam"'] 


Gl. 87-0.14 


30.2 X 27.9 


3 


50.7 


80 - 90 (90) 


G2. 11 +0.00 


31.3x26.7 


5 


79.2 


60 


Gl 1.36+0.80 


30.1 X 27.7 


3" 


48.9 


40 - 50 (40) 


G13.22-0.06 


33.1 X 32.6 


3 


76.0 


60 - 90 (80) 



Notes. The fields are named here after their approximate central Galactic coordinates (/, Z?).'*' These rms noise values refer to the maps 
used in the analysis (see text). The value given in parenthesis is the Icr noise adopted in the clumpfind analysis.''"' The observations were 
interrupted during the fourth mapping of this field, and therefore only the first three maps could be used. 



was reduced in the standard way, is clearly lower than in the 
other cases. 



2.3. C'^0(2 - 1) line observations 

From each target field, we selected seven to eight positions for 
single-pointing C'^0(2- 1) observations. These positions, which 
are listed in Table|2] were chosen from the Spitzer 8-/im images, 
and they correspond to (apparently) highly extincted parts along 
the filamentary structures near the map centres. The main pur- 
pose of making these line observations was to obtain the cloud 
radial velocity, which is needed to determine the cloud kinematic 
distance (Sect. 4.1). 

The C'^0(2 - 1) observations at 224714.199 MHz were 
carried out on 18, 22, and 26 May 2011 with APEX us- 
ing the Swedish Heterodyne Facility Instrument (SHeFI; 
Belitsky et al. 2007[ IVassilev et al. 2008al) [the heterodyne-part 
of the project 087.F-9315(A,B)]. As a frontend w e used the 
APEX-1 receiver of the SHeFI ( IVassilev et al. 2008b). The 
backend was the Fast Fourier Transfrom Spectrometer (FFTS; 
IKlein et al. 2006b with a 1 GHz bandwidth divided into 8 192 
channels. The resulting channel spacing is 122 kHz or 0.16 
km s"'. The telescope beam size (HPBW) at the observing fre- 
quency is 27'.'8. 

The observations were performed in the wobbler-switching 
mode with a 150" azimuthal throw (symmetric offsets) and a 
chopping rate of 0.5 Hz (2 s wobbler period). Total (on+off) in- 
tegration time was 5.6 min per position. The telescope pointing 
accuracy was checked by CO(2- 1) cross maps of the carbon star 
RAFGL1922, and was found to be below < 4". The focus was 
checked by measurements on Saturn. Calibration was done by 
means of the chopper-wheel technique, and the output intensity 
scale given by the system is T*^, which represents the antenna 
temperature corrected for the atmospheric attenuation. The ob- 
served intensities were converted to the main-beam brightness 
temperature scale by Tmb = T'^/jjub, where ?/mb - 0.75 is the 
main-beam efficiency at the frequency used. The single-sideband 
system temperature, in units of Tmb, was in the range 387 - 415 
K. The absolute calibration uncertainty is estimated to be around 
10%. 

The spectra were reduced using the Fortran 90 version of 
the CLASS programme from the GILD AS software package^. 
The individual spectra were averaged and the resulting spectra 
were Hanning-smoothed in order to improve the signal-to-noise 
ratio of the data. A first- or third-order polynomial was applied 
to correct the baseUne in the final spectra. The resulting 1 cr rms 



* Grenoble Image and Line Data Analysis Software is pro- 
vided and actively developed by IRAM, and is available at 



Ihttp : //www ■ iram . f r/IRAMFR/GILDAS , 



noise levels are ~ 77 - 91 mK at the smoothed resolution (4095 
channels). 

We note that the '^O nucleus has a nuclear spin of / = 5/2, 
so it has an electric quadmpole moment (-2.6 x 10"^^ cm^). The 
latter couples to the electric-field gradient at the nucleus. This 
causes the rotational lines of C'^O to have a hyperfine structure. 
The C'^0(2 - 1) line is split into nine hyperfine (hf) compo- 
nents, which cover a velocity range of about 2.36 km s"'. We 
fitted this hf structure using "method hfs" of CLASS90 to derive 
the LSR velocity (vlsr) of the emission, and FWHM linewidth 
(Av). The hf-line fitting can also be used to derive the line opti- 
cal thickness, r. However, in all spectra the hf components are 
blended together, thus the optical thickness could not be reliably 
determined. The rest frequencies and relative weights of the hf 
components were taken from Ladd et al. (1998; Table 6 therein). 

3. Observational results 

3.1. LABOCA 870-^im maps 

The obtained LABOCA maps are shown in Figs. [TJ@] In the 
right panel of each figure, we show the Spitzer 8-yum image of 
the target field, overlaid with contours of the LABOCA submm 
dust emission. 

As can be seen from the maps, the fields contain filamen- 
tary structures and clumps of different projected shapes. It can 
also be seen, especially towards the Gl .87 and G2. 1 1 fields, that 
not all 8-//m dark features are seen in submm emission. These 
may be structures with too low a column density to be detected 
with the sensitivity limit of our data. On the other hand, as was 
pointed out by Wilcock et al. (2012), some of the 8-/vm dark re- 
gions are not real dense clouds; they may just be dips in the MIR 
background that resemble the appearance of IRDCs. 

We note that in the 0.23 degr^-sized G11.36 field, the fila- 
ment in the map centre appeared to be the only submm-emitting 
object. Therefore, Fig.|3]shows only the zoomed-in view towards 
the filament. 



3.2. Clump identification 

To systematically identify the submm clumps from the 
LABOCA maps, we employed the commonly used two- 
dimensional clumpfind algorithm, clfind2d, developed by 
Williams et al. (1994). The algorithm requires two configura- 
tion parameters: i) the intensity threshold, i.e., the lowest con- 
tour level, which determines the minimum emission to be in- 
cluded into the clump; and ii) the contour level spacing, which 
determines the required "contrast" between two clumps to be 
considered as different objects. We set both of these parameters 
to the classical value of 3cr ( ~ 1 20 - 270 mJy beam" ' ), where the 
adopted 1 cr sensitivity levels are given in Col. (5) of Table[T](the 
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Table 2. Target positions of the ■'0(2 - 1) observations in the 
equatorial J2000.0 system. 



Field/ 


Q'2000.0 


^2000.0 


position 


[h:m:s] 


[":':"] 


Gl.87-0.14 






A ... 


17 50 29.6 


-27 26 02 


B ... 


17 50 31.7 


-27 25 24 


C... 


17 50 35.3 


-27 25 09 


D ... 


17 50 37.2 


-27 24 08 


E ... 


17 50 37.4 


-27 24 36 


F... 


17 50 38.0 


-27 23 02 


G... 


17 50 38.0 


-27 23 41 


G2.11+0.00 






A ... 


17 50 30.0 


-27 08 25 


B ... 


17 50 30.1 


-27 07 53 


C... 


17 50 30.3 


-27 06 37 


D ... 


17 50 30.7 


-27 07 17 


E ... 


17 50 35.6 


-27 07 15 


F... 


17 50 37.2 


-27 07 12 


G... 


17 50 38.4 


-27 07 02 


H... 


17 50 38.7 


-27 06 42 


Gll.36+0.80 






A ... 


18 07 35.0 


-18 43 51 


B ... 


18 07 35.7 


-18 42 34 


C... 


18 07 35.8 


-18 43 23 


D ... 


18 07 36.4 


-18 44 04 


E... 


18 07 36.8 


-18 41 17 


F... 


18 07 39.6 


-18 42 14 


G... 


18 07 40.5 


-18 43 16 


G13.22-0.06 






A ... 


18 14 28.2 


-17 33 28 


B ... 


18 14 31.6 


-17 32 44 


C... 


18 14 35.8 


-17 30 51 


D ... 


18 14 36.0 


-17 26 55 


E... 


18 14 36.7 


-17 29 17 


F... 


18 14 40.7 


-17 29 06 


G ... 


18 14 42.8 


-17 30 06 



value in parenthesis when the noise level varies across the map). 
Only clumps with peak flux densities greater than ~ 5cr were 
taken to be real. With these definitions, the number of clumps 
foundbyclfind2dare40, 10,7,and 34inG1.87,G2.11, 11.36, 
and G13.22, respectively. This makes the total number of identi- 
fied clumps in this survey to be 9 1 . 

In each field, the clumps are caUed SMM 1, SMM 2, etc., 
in order of increasing right ascension. The J2000.0 coordinates 
of the peak 870 yum emission, peak surface brightnesses, in- 
tegrated flux densities (within 3cr), and clump eff'ective radii 
(^eff - y/A/n, where A is the projected area within the So- 
contour) are listed in Cols. (2)-(6) of Table |3] The quoted flux 
density uncertainties are based on the rms noise values and the 
10% absolute calibration error The clump effective radii listed 
in Table[3]are not corrected for beam size. Note that the clumps 
SMM 32 and 40 in G1.87, SMM 2 in G2.11, SMM 6 in Gil. 36, 
and SMM 1, 2, and 8 in G 13. 22 are only barely resolved as their 
sizes are only slightly larger than the beam. We also note that the 
clump SMM 32 in G13.22 (Fig.|4| could be clearly resolved by 
eye into two "subclumps", but they are treated as a single source 
by cl£ind2d with our settings. 



3.3. C^'^0(2 - 1) spectra 

As shown by the green plus signs in Figs. the selected 
C'^0(2 - 1) observation target positions match very well the 



submm dust emission peaks only in G11.36. The target posi- 
tions C and D in G1.87 are within the 3o- contour of SMM 35 
(E is just outside). However, the target positions towards G2. 1 1 
do not have much, if any, associated dust emission (position 
D being an exception; it is coincident with SMM 2). Towards 
the G13.22 field the positions A, C, E, and F are rather a well 
matched with the peak positions of SMM 18, 22, 23, and 27, re- 
spectively. The obtained C'^0(2 - 1) spectra towards all target 
positions are shown in Figs. |5]-l8] Towards all fields, more than 
one velocity component is detected. This is unsurprising because 
we have observed along the Galactic midplane towards the inner 
Galaxy, where many molecular clouds along the line of sight can 
be expected. The fields G1.87, G2.11, and G13.22 show emis- 
sion at both negative and positive LSR velocities. Inspecting the 
longitude- velocity maps of CO by Dame et al. (2001), this can 
be expected at the Galactic longitudes in question. We note that 
the critical density of the C'^0(2 - 1) transition is 9.5 x 10^ 
cm"^ (assuming T = 15 K and using the data from the LAMDA 
molecular databas^ ISchoier et al. 2005T l. and therefore the line 
emission originates in dense gas. 

The C'^0(2 - 1) line parameters are given in Table |4] The 
LSR velocities and FWHM linewidths derived from hf-structure 
fits are given in Cols. (2) and (3). The peak intensities derived 
through Gaussian fitting are listed in Col. (6), and in Col. (7) 
we give the integrated line intensities computed over the veloc- 
ity range indicated in square brackets in the corresponding col- 
umn. The uncertainties in the latter two parameters take into ac- 
count the corresponding rms noise values and the 10% calibra- 
tion uncertainty. For non-detections, we provide the 3cr upper 
limit to the line intensity in Col. (6). The other parameters listed 
in Table |4] are described later. 

3.4. Clump associations 

The Spitzer 8- and 24-fj.m images were visually inspected to look 
for how the detected 870-/jm clumps appear at these MIR wave- 
lengths. The remarks concerning the 8/24 fj.m appearance of the 
clumps are given in the last column of Table |3] From the fields 
G1.87, G2.11, G11.36, and G13.22, we found, respectively, 21, 
5, 2, and 7 clumps that appear dark at both 8 and 24 jum (35 in 
total). The corresponding numbers of the clumps associated with 
either both 8 and 24 yum emission, or only with 24 yum emission, 
are 16, 3, 5, and 27 (51 in total). The type of this MIR emission 
was found to be either point sources, group of point sources, ex- 
tended, or diffuse-like. In addition, we found five clumps with 
an 8-jum point source near the submm peak position, but which 
appear dark at 24 jum. If the 8-jum source would be embedded 
within the clump, one would also expect to see emission at 24 
fim (from warm dust). As noted in Table |3] these 8-jum sources 
have5/?/fzer-GLIMPSE [3.6]-[4.5], [4.5]-[5.8], and [4.5]-[8.0] 
colours of =^ -0.03 - 0.80, 0.16 - 0.78, and 0.03 - 0.84, respec- 
tivel}0. Because these colours are not "red enough" for sources 
to be protostellar in nature, they are likely to be chance projec- 
tions of foreground stars where the emission is primarily pho- 
tospheric (cf. IGutermufli et al. 20081 IRobitaiUe et al. 2008T l. For 
sources near the Galactic plane and/or at large distances from 
the Sun, the foreground star population can be significant, and 
therefore five chance projections out of 91 clumps (5.5%) may 
not be surprising. We deal with these "foreground contaminated" 
clumps as IR-dark, although in some studies such cases are ex- 



' Ihttp : //www ■ strw . l eidenuniv . nl/~mol data/ 1 
From the GLIMPSE point source catalogue available at 
http : //irsa . ipac . caltech . edu/| 
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Jy/beam 



:G.1 87-0.1 4' 



MM 33 / 



1, C TjL 



SMM30 
SMM34 



h° SMM 27 

+ A / 



SMM 31 — 



SMM 8 



195" = 10pc 



O LABOCA beam 



SMM 7 J 
SMM 16 1^ 

1 21 _ SMM 1 



30.0 20.0 1 0.0 1 7:50:00.0 

Right ascension (J20Q0) 




30.0 20.0 1 0.0 17:50:00.0 SO.O 40:40.(> | 

Right ascension (J2000) 



Fig. 1. Left: LABOCA 870-yummap of GL87-0.14. The image is shown with hnear scaUng, and the colour bar indicates the surface- 
brightness scale in Jy beam"'. The overlaid contours go from 0.27 Jy beam"' (3cr) to 1.62 Jy beam"', in steps of 3cr. Selected clumps 
are labeled with their designation (as listed in Table|3]l. The green plus signs indicate the positions of our C'^0(2 - 1) observations 
(see Table|2]i. A scale bar indicating the 10 pc projected length is shown in the bottom left, with the assumption of a 10.57 kpc line- 
of-sight distance. The effective LABOCA beam of 19'.'9 is also shown in the lower left comer Right: Spitzer 8-//m image towards 
Gl. 87-0. 14 overlaid with LABOCA contours from the left panel. The image is shown with linear scaling, where the scale limits are 
based on the IRAF z-scale algorithm of the DS9 programme. The colour bar shows the surface-brightness scale in MJy sterad"'. 
Note that not all IR-dark regions are seen in submm emission. 




Fig. 2. Same as Fig.fflbut towards G2. 11 +0.00. The overlaid LABOCA contours go from 0.18 Jy beam"' (3o-) to 1.08 Jy beam"', 
in steps of 3cr. A scale bar indicating the 5 pc projected length is shown in the bottom left, with the assumption of a 5.51 kpc 
line-of-sight distance. The clump SMM 5, which is associated with IRAS 17474-2704, lies at a distance of 7.40 kpc (see text). 
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Fig.3. Same as Fig.[T]but towards Gil. 36+0. 80. The obtained LABOCA map is zoomed in towards the filamentary structure in 
the map centre as no submm dust emission was detected in other parts of the map. The overlaid LABOCA contours go from 0.12 
Jy beam"' (3cr) to 0.60 Jy beam"', in steps of 3cr. Note that most of the line-observation target positions are well matched with the 
submm peak positions. A scale bar indicating the 1 pc projected length is shown in the bottom left, with the assumption of a 3.27 
kpc line-of-sight distance. 



eluded from the source sample [for example, this was the case 
in the study by Chambers et al. (2009), who called such sources 
"blue cores"]. This makes the total number of IR-dark clumps in 
our survey to be 40. 

We used the SIMBAD Astronomical Databas^ to search for 
possible source associations with our clumps. The resulting as- 
sociations are given in Table |3] where we list the sources within 
about one beam size (~ 20") from the LABOCA peak position. 

In particular, 38 out of 91 clumps (42%) were found to be as- 
sociated with 1 . 1 -mm clumps from the Bolocam Galactic Plane 
Survey (BGPS; R osolowsky et al. 2010[ pLguirre et al. 201 iP * 



The BGPS survey, with an effective FWHM resolution of 33", 
has detected and catalogued about 8400 clumps. Twenty-two 
clumps (24%) are associated with SCUBA submm clumps from 
Di Francesco et al. (2008). The SCUBA maps of Di Francesco et 
al. (2008) are composed of a so-called Fundamental Map Data 
Set at 850 and 450 fim (5 061 objects), and an Extended Map 
Data Set at only 850 /urn (6 118 objects). Most of these asso- 
ciations (14 clumps) were found from the Extended Data Set 
(marked with "JCMTSE"), and eight sources have counterparts 
in the Fundamental Data Set (marked with "JCMTSF"). Thirteen 
clumps, or 14% of the sources, were found to be associated with 
IRDCs identified by Peretto & Fuller (2009E1; 

these are marked 

with "SDC" in Table [3] 

Some other associations worth mentioning here are as fol- 
lows. Four clumps were found to be associated with IRAS point 
sources: SMM 5 in G2.11 with IRAS 17474-2704, and the 
clumps SMM 8, 13, and 29 in G13.22 with IRAS 18112-1720, 



' http : // simbad . u-strasbg . fr/ simbad/ 

http : //irsa . ipac . caltech . edu/data/BOLOCAM _ GPS/| 
" The IRDC catalogue of Peretto & Fuller (2009) is available at 
www . irdarkclouds . org/ 



18114-1718, and 18117-1738, respectively. The clumps SMM 
35 and 40 in G1.87, SMM 4 in G11.36, and SMM 5, 15, and 17 
in G13.22 are associated with YSO candidates from Robitaille 
et al. (2008). Moreover, SMM 24 in G13.22 is associated with a 
candidate AGB star from Robitaille et al. (2008). Two MIR bub- 
bles from the Churchwell et al. (2006) catalogue, namely NIO 
and Nl 1, are associated with the concentration of several clumps 
in G13.22 (see Figs.|4]and[T3] and Sect. 5.1.2). Finally, we note 
that some of the clumps are associated with Class n methanol 
masers and UC Hii regions, both of which are clear signposts 
of high-mass star formation. The diffuse/extended MIR emis- 
sion seen towards some of the clumps is a typical characteristics 
of associated Hii regions and photon-dominated regions (PDRs) 
surrounding them. 

In the cases of G2.11-SMM 9, G13.22-SMM 2, G13.22- 
SMM 9, and G13.22-SMM 11, the clump appears partly asso- 
ciated (in projection) with diffuse- and/or extended-like MIR 
emission, and partly IR dark. Especially in the NlO/11 bubble 
environment, the clump classification into IR dark or IR bright 
was difficult because of the very bright and extended appearance 
of the region at 8 and 24 fim. We note that classifying clumps 
in the above four cases is a subjective process, influenced by the 
adopted colour scale of the 8- and lA-fim images. 

4. Analysis and results 

4.1. Kinematic distances 

The distance to the source is an important parameter when its 
physical properties, such as mass, are to be determined. Because 
our sources belong to the first quadrant (0° < / < 90°) in the 
inner Galaxy (i.e., interior to the solar circle), each radial velo- 
city corresponds to two kinematic-distance values along the line 
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Fig. 4. Same as Fig.[T]but towards G13.22-0.06. The overlaid LABOCA contours go from 0.24 Jy beam"' (3cr) to 1.44 Jy beam"', 
in steps of 3cr. The field contains several clumps at different kinematic distances, and therefore only the angular scale bar is shown. 
The green circles in the Spitzer 8 fim image on the right panel indicate the positions and outer radii of the bubbles NIO and Nil 
from Churchwell et al. (2006; their Table 2). 



of sight. However, in the first Galactic quadrant, the radial ve- 
locity of the source increases as a function of distance up to the 
tangent point. At this point, the source's velocity vector and the 
line of sight are aligned with each other, and the radial velocity 
has its maximum value. After passing the tangent point, the ra- 
dial velocity starts to decrease as a function of distance, all the 
way down to negative values dRoman-Duval et al. 2 009; Fig. 2 
therein). In general, for sources associated with IRDCs, it is rea- 
sonable to assume, and often has been assumed, that they lie at 
the near distance, because in that case there is more background 
radiation against which to see the cloud in absorption. 

To calculate the kinematic distances of G1.87, G2.11, and 
Gil. 36, we adopted the average C'^0(2- 1) radial velocities to- 
wards each field (using the velocities at which most lines were 
seen, and excluding the additional velocity components). The 
obtained average LSR velocities for G1.87, G2.11, and G11.36 
are, respectively, -41 .5, 16.7 (from positions A and B), and 27.8 
km s-i. The clump SMM 5 (IRAS 17474-2704) in G2.11 is 
an exception. Although being near to our line observation po- 
sitions in the plane of the sky, it has the peak Class ii CH3OH 
maser emission at 63 km s~' jCaswell et al. 19931 l. The OH and 
H2O masers towards this source also peak at comparable veloci- 
ties dForster & Caswell 1989l l. Therefore, for G2.1 1-SMM 5 we 
adopted the velocity 63 km s ' . 

The average C'^0(2 - 1) LSR velocity for G13.22 is 35.9 
km s ' . However, from the HCO^ and N2H^ survey of the BGPS 
1.1 -mm clumps by Schlingman et al. (2011), we obtained velo- 
city information for the following clumps in G13.22: SMM 5, 
10, 13, 17, 24, 25, 27, 28, 29, and 32. As shown in Fig. S the 
other clumps in G13.22 are seen in projection close to the above 
listed clumps. For example, the clumps near the bubbles NlO/1 1, 
such as SMM 5 and 7, are likely physically connected, and SMM 
31 is likely associated with SMM 28. Another example is SMM 



26, which is likely a member of the filament connecting SMM 
24 and 29. Therefore, for the field G13.22 we do not adopt our 
C'^O-derived radial velocity, but instead use the values from 
Schlingman et al. (201 1). Sewilo et al. (2004), using the HI IQa 
and H2CO observations, were able to distinguish between the 
neai- and far distances of G13.22-SMM 29 (IRAS 181 17-1738); 
H2CO absorption was seen between the source velocity and the 
velocity at the tangent point, placing the source at the far dis- 
tances 

We employed the rotation curve of the Galaxy by Reid et 
al. (2009), which is based on direct measurements of trigono- 
metric parallaxes and the proper motions of masers in high-mass 
star-forming regions. The best-fit rotation parameters of Reid et 
al. (2009) are (0o, ^o)=(254 km s"', 8.4 kpc), where 0o is the 
orbital velocity of the Sun around the Galactic Centre, and Rq 
is the solar galactocentric distance. The resulting near and far 
kinematic distances, Jneai and lifai, and the Galactocentric dis- 
tances, Rcc, are given in Table |5] Unless otherwise stated, the 
near distance was adopted. We note that Rqc does not have a 
distance ambiguity. All the other clumps in G13.22 for which we 
adopted the distances derived using the data from Schlingman et 
al. (201 1) are listed in Col. (5) of Table |5] We note that our dis- 
tances for the G13.22 clumps correspond to the values reported 
by Schlingman et al. (201 1) who also used the Reid et al. (2009) 
rotation curve. 

As discussed above, the negative radial velocity of G1.87 
places it at the far distance. In principle, it could also belong 
to the near 3 -kpc arm at d ~ 5 kpc dPame & Thaddeus 20081 
IGreen et al. 201 II Tackenberget al. 2012 , but the far solution 



is adopted in this work. Another noteworthy issue to raise is 
the fact that our line observations were only made towards the 



We note that the HllOa velocity of 40.6 ± 3.8 km s ' is comparable 
to the HCO+ velocity of 44.64 km s"' from Schlingman et al. (2011). 
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Fig. 5. Smoothed C'^0(2 - 1) spectra towards selected positions in Gl. 87-0. 14. Most lines are seen at 41 km s '. Hyperfine- 

structure fits to the lines are overlaid in green. No line emission is detected towards position B, whereas two velocity components 
are seen towards positions C and F. 



filamentary structures near the map centres, and therefore the 
derived radial velocities may not apply for all the clumps de- 
tected in the field. This may especially be the case towards 
G1.87, where many clumps are detected at low Galactic longi- 
tudes. Additional distance-uncertainty towards G1.87 is caused 
by the possible association of some of the clumps with the near 
3-kpc arm. On the other hand, the field Gil. 36, where only one 
velocity-coherent filament is detected, can be considered to have 
the most reliable distance estimate among our target fields. 



4.2. Temperatures 

The dust and gas temperatures, Tdust and Tkin, are also essential 
knowledge when studying the physics and chemistry of molecu- 
lar clumps. To our knowledge, for only one clump in our sample, 
namely SMM 5 in G13.22, the gas kinetic temperature measure- 
ment has been published. Pillai et al. (2007) derived the NH3 ro- 
tation temperature of Ti-ot - 17.4± 1 .3 K for SMM 5 (their source 
G13. 18+0.06). Using the T^oi - r^n relationship from Tafalla et 
al. (2004), we obtain Tun = 20.5 + 1.7 kQ. We note that at 
high densities of n(H2) > 3 x 10"* cm"^, where collisional cou- 
pling between the gas and dust becomes efficient, the gas and 
dust temperatures are expected to be similar, T[^„ ^ Tdust (e.g., 
IGaUi et al. 20"02l) . 

For those four clumps in our sample which are associated 
with IRAS point sources, we estimated the dust temperature to 
be the same as the 60/100-yum colour temperature defined by 
Henning et al. (1990) as 



" Pillai et al. (2007) only reported the rrot(NH3) values. We note that 
the NH3 radial velocity measured by Pillai et al. (2007) is in good agree- 
ment with the HCO^ radial velocity from Schlingman et al. (201 1). 



rdus.-re|^|-96 



(1) 



In this formula, /3 is the dust emissivity index, and 5^ is the flux 
density at the wavelength A. The value of /3 was set to be 1.8 to 
be consistent with the Ossenkopf & Henning (1994) dust model 
discussed in Sect. 4.3. For the IRAS sources 17474-2704, 18112- 
1720, 18114-1718, and 18117-1738, i.e., for the clumps G2.11- 
SMM 5, G13.22-SMM 8, G13.22-SMM 13, and G13.22-SMM 
29, we derived the Tdust values of 30.0, 18.9, 35.5, and 18.9 K, 
respectively. 

For the remaining 86 clumps the dust temperatures we 
assumed to be the following: 15 K for IR-dark clumps, 20 
K for IR-bright clumps (8/24 fim emission), and 30 K for 
clumps associated with Hii regions/radio sources. The choice 
of Tdust = 15 K or 20 K for most of our clumps is expected to 
be reasonable, because previous molecular-line observations of 
clumps within IRDCs have shown the typical gas kinetic tem- 
perature to lie in the range r^n ~ 10 - 20 K ( Carey et al. 1998[ 
Teyssier et al. 2002| [Sridhara n eLaL20 05J IPillai et al."2006t 
Sakaietal. 2008: Zhang et al. 201 1[ IDevine et al. 201 U 



Ragan et al. 201 1) . For comparison, in their study of a massive 
clump associated with an IRDC, Hennemann et al. (2009) 
derived the Tdust values of 22 K and 15 K for the clump's 
substructures ('cores') with and without 24-jum sources, respec- 
tively. Rathbome et al. (2010) also found that 24-fim bright 
clumps embedded in IRDCs are warmer than their 24-fim dark 
counterparts. The assumption of a slightly higher temperature of 
30 K in clumps with embedded Hii regions is supported by the 
30-K temperature of the UC Hii region G2.11-SMM 5 derived 
above, and also consistent with some other observational results 



(see, e.g., Sreenilayam & Fich 201 1 and references therein) 
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Fig. 6. Same as Fig.|5]but towards the selected positions in G2.1 1+0.00. Most detected lines are near ~ 16 km s No lines were 
detected towards positions D, F, and G. Two velocity components are seen towards position C. 



Finally, we note that the BGPS 1.1 -mm data available for 
42% of our clumps could, in principle, be used to estimate the 
dust colour temperature. However, there are three factors that 
would hamper this analysis: i) the Bolocam 1.1 -mm flux den- 
sities may be somewhat uncertain (due to spatial filtering and 
caUbration issues), and should be multiplied by 1.5 ± 0.15 as 
recommended by Aguirre et al. (201 1); ii) one should make as- 
sumptions about the value of /? as above; and Hi) the LABOCA 
and Bolocam wavelengths are quite close to each other (0.87 mm 
vs. 1.1 mm), and therefore the corresponding flux densities are 
comparable to each other In addition, our LABOCA data should 
be smoothed to the 33" BGPS resolution for a proper compari- 
son. For these reasons, it seems justified to assume that the Tdust 
values are similar to those observed in some other sources of 
similar type. 

4.3. Radii, masses, and H2 column and number densities 

The linear clump effective radii (in pc) were derived from the 
angular radii and kinematic distances as ^efF[pc] = 7?eff [rad] x 
c/[pc]. The error in R^ff was computed from the average value 
between the +-distance errors quoted in Table |5] 

The clump masses, M, over an effective area of radius 
Rgff were estimated from the integrated 870-/zm flux densities 
and using the standard optically thin dust emission formulation 
[see, e.g., Eq. (6) in [Miettinen & Harju 20T0t hereafter MH 10]. 



Following Hatchell et al. (2007), the peak surface brightness 
was used if it was higher than the integrated flux density (this 
was the case for G1.87-SMM 40, G2.1 1-SMM 2, Gl 1.36-SMM 
6, and G13.22-SMM 8; all but G2.1 1-SMM 2 are associated 
with previously known sources, and are therefore likely to be 
real). The distances and dust temperatures used were as ex- 
plained above. The dust opacity per unit dust mass at 870 fim 
was taken to be /cgvo - 1.38 cm^ g"'. This value was extrapo- 



lated from the Ossenkopf & Henning (1994) model describing 
graphite-silicate dust grains that have coagulated and accreted 
thin ice mantles over a period of 10^ yr at a gas density of 
riH = n(H) + 2n(H2) ^ 2n(H2) = 10^ cm^^Qln this model, the 
dust emissivity index is yS ^ 1.8, as determined from the slope 
between 350 and 1300 fim (k^ cc A^^). For the average dust- 
to-gas mass ratio, Ra = (Mdust/A^gas), we adopted the canonical 
value 1/100 (Spitzer 1978jl. The uncertainty in mass was propa- 
gated from the uncertainties in flux density and source distance 
(the 1.7-K temperature error was also employed in the case of 
G13.22-SMM5). 

The peak beam-averaged H2 column densities, A^(H2), were 
computed from the peak surface brightnesses in a standard way 
[see, e.g., Eq. (8) in MHIO]. For this calculation we assumed a 
He/H abundance ratio of 0.1, which leads to the mean molecular 
weight per H2 molecule of yUe, - 2.8. Other parameters (Tdust, 
Km), Rd) were the same as used in the mass calculation. The error 
in A^(H2) is solely based on the uncertainty in the peak surface 
brightness (except in the case of G13.22-SMM 5, where Tdust 
error was also used). The volume-averaged H2 number densities 
over 7?efF, (n(H2)), were calculated using Eq. (7) of MHIO, and 
the corresponding errors were propagated from those of M and 

The values of the physical parameters derived above are 
listed in Cols. (6)-(9) of Table|3] Their distributions are shown in 
Fig-El panels a)-d) show the histograms of 7?eif, M, N(li2), and 
(n(H2)), respectively. The parameters' mean and median values 
are described in the figure caption. 



" We note that in the ATLASGAL study of Schuller et al. (2009), the 
value Ksjo = 1.85 cm^ g"' was used. 
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Fig. 7. Same as Fig.|5]but towards the selected positions in Gil. 36+0. 80. We note that clear line emission at ~ 28 km s"' is seen 
towards all positions. Two velocity components are detected towards positions A, B and F, but the secondary line is very weak in 
the former two case. 



4.4. Non-thermal velocity dispersion 

The measured C'^0(2 - 1) linewidths were used to calculate 
the non-thermal portion of the line-of-sight velocity dispersion, 
ctnt (averaged over a 27'.'8 beam); see, e.g., Eq. (1) in Miettinen 
(2012, hereafter M12). We also computed the ratio ctnt/cs, 
where Cs = -y/A:B7'kin/jUp'WH is the isothermal sound speed with 
the Boltzmann constant, fip - 2.33 the mean molecular weight 
per free particle for solar composition (He/H- 0.1), and niu 
the mass of the hydrogen atom. Some of our line observations 
are probing embedded YSOs, and = 20 K was assumed 
for these accordingly. For other positions, the value = 15 
K was adopted. The values of crm and crt^j/cs are listed in 
Cols. (4) and (5) of Table |4] The uncertainties are based on 
the linewidth uncertainties. The vast majority of the positions 
show supersonic non-thermal motions (ctnx > Cs), which are 
presumably due to turbulence. This is a general observational 
feature of the interstellar molecular clouds (e.g., Larson 19811 
iMcKee & Zweibel 1992: Heyer & Brunt 2004[ ). 

4.5. column densities, fractional abundances, and CO 
depletion factors 

The beam-averaged column densities, A^(C'^O), were de- 
rived following the standard LTE analysis outlined, e.g., in the 
paper by M12 (Appendix A. 3 therein). In brief, we have assumed 
optically thin line emission, and computed the A^(C'^O) values 
from the integrated line intensities. In the two cases where the 
detected line does not cover all the hf components (see Table |4|, 
the column densities were scaled by the inverse of the rela- 
tive line strength within the detected line. The C'^0(2 - 1) 
transition was assumed to be thermalised at the gas tempe- 
rature of the target position as seen in previous studies (e.g., 
IMiettinen et al. 201 11 1. Therefore, the line excitation temperature 



was adopted to be Tex = 15 K towards all the positions except 
those associated with MIR emission, where Tex = 20 K was 
used. We note that if Tex increases from 10 to 20 K, the column 
density decreases by a factor of 2.3 (or by a factor of 1 .3 for Tgx 
of 15-20 K). The errors in A^(C'^O) were propagated from those 
associated with the integrated intensity. 

The fractional C'^O abundances were computed by divid- 
ing the C'^O column density by the H2 column density as 
ji:(Ci^O) = NiC"0)/Nin2). For this puipose, the NiHjJ values 
were derived from the LABOCA maps smoothed to the 27'.' 8 
resolution of the line observations. The uncertainty in x(C'^O) 
was derived from both the errors in the C'^O and H2 column 
densities. 

The CO depletion factors, /d, were also derived following 
the analysis in M12 (see Appendix A. 4 and references therein). 
In summary, the Galactocentric distance of the source was first 
used to estimate the appropriate canonical (or undepleted) CO 
abundance. For example, the "field" Rcc values shown in Table|5] 
lead to the CO abundances of x(CO) - 1.5 - 2.2 x 10 
We also employd the /?gc -dependent ['^0]/['^0] ratio, which 
is ^ 166 - 345 for the "field" Rcc values, and adopted the 
["^0]/['^0] ratio of 3.52. This way, we estimated the canoni- 
cal C'^O abundance and calculated the depletion factor as /d = 
x(C'^0)can/.^(C'^0)obs- The /d uncertainty was propagated from 
that in the observed fractional abundance. The results of the cal- 
culations presented in this section are shown in Table |4](the last 
three columns). We note that we only report the values of /d 
for the "main" velocity components, because the line profiles in 
these cases are more reliable than those of the additional velocity 
components. We also stress that the formal /d errors probably 
underestimate the true uncertainties by a factor of > 2 - 3 be- 
cause of the uncertainties in the assumptions used (Tex, oxygen- 
isotopic ratios, etc.). Therefore, some of the /d values might be 
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Fig. 8. Same as Fig.|5]but towards the selected positions in G13.22-0.06. We note that clear line emission at ~ 37 km s ' is seen 
towards aU positions. Two velocity components are detected towards positions A-C and G. 



significantly lower than reported here. In this regard, the appar- 
ent variation of /d between different sources might not be robust. 

4.6. Analysis of the filamentary IRDC G1 1.36+0.80 

The G11.36 cloud represents the best example of a filamentary 
IRDC in our survey. In addition, as demonstrated in Fig. |3] our 
C'^O line observations probe the clumps along the filament. For 
these reasons, we will analyse the filament's properties in more 
detail here. 

The north-south oriented G11.36 filament has a total pro- 
jected length of about 4.05' or 3.85 pc, and its average radius 
is about 0.4 pc. By excluding the clumps SMM 6 and 7, which 
lie on the side of the filament, we can estimate the mass of the 
filament to be ~ 762 M0 as the sum of the clump masses within 
it. This makes the filament's mass per unit length, or line mass, 
to be Miine ~ 198 M0 pc"^ The average projected separation 
between the five clumps in the filament is about 0.9 pc. 

As shown in Col. (5) of Table [H the filament appears to 
be dominated by supersonic non-thermal motions (o-t^j/c^ ^ 
1.5 - 3.7). Therefore, to examine the dynamical state of the fi- 
lament, we calculate its virial mass per unit length as M™^ - 
2{cr^}/G, where (cr^) is the square of the total (thermal-i-non- 
thermal) velocity dispersion, i.e., the square of the effec- 
tive sound speed (Ceff), and G is the gravitational constant 
(Fiege & Pudritz 2000aV Assuming that the gas kinetic temper- 
ature is Tkin = 15 K, and using the average non-thermal velocity 
dispersion of 0.69 km s ', we obtain Ceff = 0.73 km s"^ Thus, 
we derive the values M^^^ ^ 246 M© pc"' and Miine/M™^ ^ 0.8. 
This impUes that Gl 1 .36 as a whole is close to virial equilibrium. 

The magnetohydrodynamic "sausage"-type instability the- 
ory predicts that growing perturbations can fragment a self- 
gravitating fluid cylinder into successive condensations with al- 
most periodic separations (e.g., IChandrasekhar & Fermi 19531 



[Nagasawa 1987[ see also IJackson et al. 201 Oi l. This separation 
distance corresponds to the wavelength of the fastest growing 
mode, /Imax - (^max being the wavenumber), which, 

in turn, is twice the wavelength of the axisymmetric perturba- 
tions to which the cylinder is unstable. Jackson et al. (2010) 
and Ml 2 applied the sausage-type instability predictions about 
clump separations to the filamentary IRDCs "Nessie" Nebula 
and G304. 74-1-01. 32, respectively. Both studies found that the 
observed clump separations agree with theoretical predictions 
when non-thermal (turbulent) motions are taken into account. In 
this case, the fastest growing mode in an isothermal, infinitely 
long gas cylinder appears at 



22//, 



eif 



22 X 



Ceff 



V47rGpo 



6.2 X 



Ceff 



(2) 



where H^a is the effective radial scale height with po the cen- 
tral gas-mass density along the cylinder's axis. If we compute 
Po assuming that the central number density is 10^ cm""*, which 
seems reasonable because even the volume-average densities are 
~ 10"* cm""* (Table O, we derive the values //eff ^ 0.04 pc and 
^max - 0.9 pc. This is in excellent agreement with the observed 
average clump separation of ~ 0.9 pc. The cloud may therefore 
lie close to the plane of the sky. Moreover, the wavelength of the 
fastest growing perturbation appears to be approximately equal 
to the filament's diameter of ~ 0.8 pc, in accordance with theory 
(INakamura et al. 1993l l. 

In this theoretical framework, the clump masses should be 
less than M ~ /ImaxMine (see IJackson et al. 20101 l. Using either 
the observed line mass or the MJ^^^ value, the predicted maxi- 
mum mass is about 180 or 220 Mq, respectively. The estimated 
clump masses in the filament are in the range ~ 101 - 192 M0, 
which conforms to the theoretical expectation. The fragmenta- 
tion timescale for a filament of radius Rfn is expected to be com- 
parable to its radial signal crossing time, Tcross = ^fii/(o"^)'^^ 
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Table 4. C 0(2 - 1) line parameters, column densities and fractional abundances, and CO depletion factors. 



Position^' 


vlsr 


Av 






7'mb 


J TMBdv* 


N(C'''0) 


x(C"0) 


/d(CO) 




[kms-'] 


[kms-'] 


[kms-'] 


[K] 


[Kkms-'] 


[10" cm-=] 


[10-8] 




Gl.87-0.14 




















A... 


-40.5 ±0.2 


1.36 ±0.58 


0.57 ±0.25 


2.5 ± 1.1 


0.24 ± 0.06 


0.63 ±0.13 [-42.81. J7.90] 


3.4 ± 0.7 


7.4 ±4.0 


5.0 ±2.7 


B ... 










< 0.26 










C ... 


-41.2 ±0.4 


3.98 ± 1.11 


1.69 ±0.47 


6.3 ± 1.8 


0.20 ± 0.06 


0.81 ±0.19 [-45.15, -38.14] 


2.9 ± 0.7 


3.1 ± 1.0 


11.9±3.8 


C(2nd) . . . 


47.0 ± 0.5 


2.23 ± 0.61 


0.94 ± 0.26 


3.5 ± 1.0 


0.14 ±0.06 


0.70 ±0.18 [42.58, 49.60] 


2.7 ±0.5 


2.9 ±0.8 




D ... 


-4I.4±0.I 


0.55 ± 0.44 


0.22 ± 0.20 


0.8 ± 0.7 


0.34 ± 0.07 


0.65 ±0.11 [-43.98,-40.01] 


2.8 ±0.5 


2.1 ±0.5 


17.6 ±4.2 


E... 


-41.6 ±0.2 


1.54 ±0.45 


0.65 ±0.19 


2.8 ±0.8 


0.30 ± 0.07 


0.68 ±0.12 [-44.22, -39.54] 


3.7 ± 0.6 


2.1 ±0.5 


17.6 ±4.2 


F... 


-41.3 ±0.2 


0.81 ±0.76 


0.34 ±0.33 


1.5 ± 1.4 


0.31 ±0.05 


0.43 ±0.10 [-44.22, -40.94] 


2.3 ±0.5 


1.9 ±0.6 


19.4 ±6.1 


F(2nd) . . . 


61.5 ±0.3 


1.59 ±0.84 


0.67 ±0.36 


2.9 ± 1.6 


0.18 ±0.05 


0.39 ±0.11 [59.19, 63.17] 


2.1 ±0.6 


1.8 ±0.6 




G... 


-41.9 ±0.1 


1.39 ±0.62 


0.59 ±0.26 


2.5 ± 1.1 


0.28 ± 0.08 


0.75 ±0.18 [-45.39, -39.07] 


4.1 ± 1.0 


4.6 ± 1.7 


8.0 ±3.0 


G2.11+0.00 




















A... 


16.8 ±0.1 


0.83 ± 0.24 


0.35 ±0.10 


1.5 ±0.5 


0.35 ± 0.04 


0.37 ±0.08 [15.82, 18.16] 


2.0 ± 0.5 


4.9 ± 2.2 


5.5 ±2.5 


B ... 


16.5 ±0.3 


1.33 ±0.47 


0.56 ±0.20 


2.4 ±0.9 


0.19 ±0.07 


0.28 ±0.08 [14.88, 18.16] 


1.5 ±0.5 


7.0 ±5.5 


3.8 ±3.0 


C ... 


15.7 ±0.3 


1.43 ±0.40 


0.60 ±0.17 


2.6 ±0.7 


0.10 ±0.06 


0.37 ±0.13 [13.13, 17.92] 


2.0 ±0.7 


4.4 ±2.1 


6.1 ±2.9 


C(2nd) . . . 


-11.5±0.1 


0.57 ±0.41 


0.23 ±0.18 


1.0 ±0.8 


0.25 ± 0.04 


0.22 ± 0.07 [-12.43, -10.43; 95.1%] 


1.2 ±0.4 


2.6 ± 1.2 




D ... 










<0.26 










E... 


-19.4 ±0.3 


2.04 ± 0.76 


0.86 ± 0.32 


3.7 ± 1.4 


0.16 ±0.07 


0.37 ±0.10 [-20.91, -17.52] 


2.0 ±0.6 


6.2 ±3.4 


4.3 ± 2.4 


F... 










<0.27 










G... 










<0.23 










H ... 


23.8 ± 0.2 


0.85 ± 0.42 


0.36 ±0.18 


1.5 ±0.8 


0.17 ±0.05 


0.25 ± 0.09 [22.02, 25.53] 


1.4 ±0.5 


1.3 ±0.5 


20.7 ± 8.0 


Gll.36+0.80 




















A... 


27.6 ±0.1 


1.93 ±0.34 


0.82 ±0.14 


3.5 ±0.6 


0.42 ± 0.09 


1.05 ±0.12 [24.39, 30.47] 


4.5 ± 0.6 


3.1 ±0.5 


3.9 ±0.6 


A(2nd) . . . 


56.4 ± 0.2 


0.62 ±0.17 


0.25 ± 0.08 


0.9 ± 0.3 


0.16 ±0.05 


0.31 ±0.12 [54.25, 58.31] 


1.7 ±0.6 


1.2 ±0.4 




B ... 


28.0 ±0.1 


1.06 ±0.12 


0.45 ± 0.05 


1.9 ±0.2 


0.88 ±0.16 


2.04 ± 0.23 [25.09, 30.70] 


11.0± 1.2 


7.3 ± 1.1 


1.6 ±0.2 


B(2nd) . . . 


55.8 ±0.2 


0.81 ±0.59 


0.34 ±0.26 


1.5 ± 1.1 


0.18 ±0.04 


0.16 ±0.10 [54.33, 57.37] 


0.9 ± 0.5 


0.6 ±0.3 




C... 


27.7 ±0.1 


2.01 ±0.21 


0.85 ± 0.09 


3.7 ± 0.4 


0.68 ±0.11 


1.73 ±0.21 [23.73, 30.41] 


9.3 ± 1.1 


4.8 ±0.8 


2.5 ± 0.4 


D ... 


27.8 ±0.1 


1.95 ±0.21 


0.83 ± 0.09 


3.6 ±0.4 


0.62 ±0.11 


1.55 ±0.19 [24.62, 30.70] 


8.4 ± 1.0 


5.6 ±0.9 


2.1 ±0.3 


E... 


27.1 ±0.1 


1.18 ±0.25 


0.50 ±0.11 


2.2 ±0.5 


0.53 ± 0.09 


1.11 ±0.15 [24.39, 29.53] 


5.4 ±0.7 


3.8 ±0.7 


3.2 ±0.6 


F... 


27.9 ±0.1 


1.50 ±0.45 


0.63 ±0.19 


2.7 ±0.8 


0.59 ± 0.07 


1.32 ±0.17 [25.56, 30.47] 


7.1 ±0.9 


12.5 ±3.2 


I.0±0.2 


F(2nd) . . . 


83.7 ±0.1 


0.54 ±0.13 


0.22 ± 0.06 


I.0±0.2 


0.43 ± 0.06 


0.30 ± 0.07 [82.64, 84.98; 84.3%] 


1.9 ±0.5 


3.3 ± 1.1 




G... 


28.6 ±0.1 


0.71 ±0.09 


0.29 ± 0.04 


1.3 ±0.2 


0.64 ±0.15 


1.47 ±0.19 [26.02, 31.17] 


6.0 ±0.8 


8.2 ± 1.5 


1.5 ±0.3 


G13.22-0.06 




















A... 


37.0 ±0.1 


1.06 ±0.17 


0.44 ± 0.07 


1.7 ±0.3 


0.47 ± 0.06 


0.69 ±0.11 [34.45, 38.66] 


2.8 ± 0.5 


1.4 ±0.3 


12.0 ±2.6 


A(2nd) . . . 


15.1 ±0.1 


0.69 ±0.16 


0.29 ± 0.07 


1.2 ±0.3 


0.35 ± 0.05 


0.47 ±0.08 [12.45, 16.90] 


1.9 ±0.3 


0.9 ± 0.2 




B ... 


35.9 ± 0.3 


1.11 ± 1.08 


0.47 ± 0.46 


2.0 ±2.0 


0.33 ± 0.08 


0.56 ±0.10 [33.74, 37.02] 


3.0 ±0.5 


1.5 ±0.4 


8.5 ±2.3 


B(2nd) . . . 


54.0 ± 0.3 


1.05 ± 1.08 


0.44 ± 0.46 


1.9 ±2.0 


0.48 ± 0.09 


1.02 ±0.15 [51.29, 56.91] 


5.5 ±0.8 


2.8 ±0.6 




C... 


36.6 ±0.1 


1.71 ±0.30 


0.72 ±0.13 


3.1 ±0.6 


0.46 ±0.10 


1.07 ±0.16 [33.51, 38.89] 


5.8 ±0.9 


3.4 ±0.8 


3.8 ± 0.9 


C(2nd) . . . 


54.9 ± 0.2 


1.12±0.19 


0.47 ± 0.08 


2.0 ±0.4 


0.27 ± 0.09 


0.77 ±0.14 [51.99, 56.67] 


4.2 ±0.8 


2.5 ±0.6 




D... 


37.3 ±0.1 


1. 17 ±0.11 


0.49 ± 0.05 


2.1 ±0.2 


0.92 ±0.12 


1.69 ±0.20 [34.68, 39.59] 


9.1 ± 1.1 


12.3 ±4.1 


I.0±0.3 


E... 


38.0 ±0.2 


2.58 ±0.66 


1.09 ±0.28 


4.1 ± 1.1 


0.43 ± 0.08 


1.41 ±0.22 [34.91,41.46] 


5.8 ±0.9 


3.9 ±0.8 


3.3 ±0.7 


F... 


36.7 ±0.1 


1.47 ±0.29 


0.62 ±0.12 


2.3 ± 0.5 


0.62 ±0.14 


1.88 ±0.23 [33.51,40.06] 


7.7 ±0.9 


1.3 ±0.2 


9.9 ± 1.5 


G... 


35.9 ±0.3 


0.87 ± 1.08 


0.36 ± 0.47 


1.6 ±2.0 


0.46 ± 0.07 


0.75 ±0.12 [33.98, 37.96] 


4.1 ±0.6 


4.0 ± 1.2 


3.2 ± I.O 


G(2nd) . . . 


-5.7 ± 0.6 


2.06 ± 1.31 


0.87 ±0.56 


3.8 ±2.4 


0.16 ±0.06 


0.57 ±0.15 [-9.20, -3.11] 


3.1 ±0.8 


3.1 ± 1.1 





Notes. Columns (2)-(10) of this table are as follows: (2) LSR velocity; (3) FWHM linewidth; (4) one-dimensional non-thermal velocity 
dispersion; (5) the ratio of the non-thermal velocity dispersion to the isothermal sound speed; (6) peak line intensity; (7) integrated intensity; 
(8) total column density; (9) fractional abundance of C"0; (10) CO depletion factor (only towards the 'main' velocity components).*"' The 
additional velocity components are labeled with "(2nd)".'*' The velocity range used in the calculation is given in square brackets. The percentage 
value for G2. 1 1-C(2nd) and Gl 1 .36-F(2nd) indicates the contribution of hf components' intensity lying within the detected line. 



[see Eq. (26) in |Fiege & Pudritz 2000b| . For Gil. 36, this is es- 
timated to be ~ 5.4 X 10^ yr. 



4. 7. Virial analysis of the clumps in G11 .36+0.80 

We employed the C"0(2 - 1) Hnewidths and Eq. (3) in M12 
to calculate the virial masses, Mvir, of the clumps in G11.36. It 
was assumed that the clumps have a density profile of the form 
n(r) oc r"' " found by, e.g., Beuther et al. (2002) for high-mass 
star-forming clumps. The corresponding virial parameters were 
derived following the definition of Bertoldi & McKee (1992), 
i.e., ttvir = Mvii /M. If ffvir = 1, the clump is in virial equilibrium, 
whereas the clumps with avii ^ 2 are taken to gravitationally 
bound. The deiived values of Mvi, and ffvir are given in Table |6] 
and ffvir is also plotted as a function of the clump mass in Fig.fTOl 
The error in Mvir was propagated from those associated with Av 
and Rgff, whereas the error of ayir includes the uncertainties in 
both the mass values Mvir and M. Although the unceitainties 
are large, it seems that five out of seven clumps in G11.36 are 
gravitationally bound. The coiTesponding ratio is 4/5 when con- 
sidering only the clumps in the filament We note that the submm 
peak of SMM 4 was missed by the line observations, position 



Table 6. Virial masses and parameters for the clumps in Gl 1 .36. 



Source" 








[Mq] 





SMM 1(A) 


253 


+ 90 


2.0 ± 


0.9 


SMM 2(C) 


259 


±66 


1.7 ± 


0.7 


SMM 3(B) 


115 


±26 


0.6 ± 


0.2 


SMM 4(D) 


258 


±64 


2.6 ± 


1.0 


SMM 5(E) 


150 


±55 


0.8 ± 


0.4 


SMM 6(F) 


74: 


t41 


3.0 ± 


2.0 


SMM 7(G) 


49: 


t 10 


0.9 ± 


0.3 



Notes. '"' The letter in parenthesis indicates the line observation position 
as shown in Fig. [3] Note that position D is not very well coincident with 
the submm peak of SMM 4. 



D being closest to it. Therefore, the virial parameter of SMM 4 
should be interpreted with some caution. 
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Field/ 






Rgc 


Remark^ 


clump* 


[kpc] 

,^ w 


[kpc] 

in c-7+0.5«c 


[kpc] 





G2.11+0.00 

SMM5 
Gll.36+0.80 
G13.22-0.06 

SMM5 
SMM 10 
SMM 13 
SMM 17 
SMM 24 
SMM 25 
SMM 27 
SMM 28 
SMM 29 
SMM 32 



5.51+"-" 
7 40+^:1^ 

'■^"-0.12 

2 27+0-47 

3.54+''-^^ 
4 24+8-^^ 
4 

^■^^-0.34 

(3.98^11) 

3 56+^-^^ 
■'■J°-Q.44 

446+B.23 
(4.03+li) 



Q ^Q+0.12 

'■•"-0.10 
r-1+0.56 



4.4r 



37 



13.20 
12.82_ 
12.12+ 
12.04+ 
14.39+ 
12.80+ 
12.38+5; 
11.88+' 
12.79+ 
11.89+ 
12.33+5; 
11.94 ' 



-047 
+045 



2.90 
1.04 
5.23 
5.02 
4.38 
4.32 
6.52 
5.00 
4.61 
4.17 
5.00 
4.19 
4.57 
4.23 



- for all clumps except SMM 5 



this dnear was not used for any of the clumps 
rfnear also for SMM 2, 4, 6, 7, 9, 11, 12 
rf„Mr also for SMM 1,3,8 

d„Mr also for SMM 15, 20 

also for SMM 26, 33, 34 
rfnear also for SMM 18,21 
rfnear also for SMM 22, 23 
rf„ear also for SMM 14, 16, 19, 30, 31 



Notes. The near and far Idnematic distances were determined using the Reid et al. (2009) rotation curve of the Galaxy. The near distance 
was adopted unless otherwise stated.**' The fields G2.11 and G13.22 are known to contain sources at different distances than determined here 
from our radial velocities (see text).'^' Other clumps for which the quoted distance was adopted.'''' The negative LSR velocity of G1.87 
places it at the far distance. However, it could also be a member of the near 3-kpc arm. The near distance derived here is comparable to the ~ 5 
kpc distance to the near 3-kpc arm.''"' Far distance was adopted. 



5 - 



4 r 



■5 ^ r 
a' ; 




f , , , , I , , , , I , , , , I , , , , I , , , , I ^ 
50 100 150 200 250 

M [M 1 

Fig. 10. Virial parameter vs. mass for the clumps in Gl 1 .36. The 
dashed line indicates the virial-equilibrium limit of Ofvir = 1, and 
the dash-dotted line shows the limit of gravitational boundedness 

or ffvir - 2. 



5. Discussion 

5. 7. On the nature of the detected clumps 

The sizes and masses of the clumps detected in this study indi- 
cate that they are the precursors/formation sites of stellar clusters 
and groups, rather than those of single stars or low-order multi- 
ples. Among our sample, 40 clumps are IR dark, and 51 clumps 
are IR bright, leading to the relative fractions of 44% and 56%, 
respectively. Some of the detected IR-dark clumps may represent 
the high-mass starless "cores" (HMSCs). However, it is fully 
possible that at least some of them do contain embedded YSOs, 
but which are not bright enough to be detected with the present 
sensitivity. The clumps associated with point-like MIR emission 
at 8 and 24 yum, or only at 24 //m, are likely to host YSOs. The 
only exceptions may be the cases where the clump is partly IR 
dark, and partly associated with extended MIR emission, possi- 
bly from the nearby sources as in the case of the NlO/11 bub- 
ble environment. For example, star-formation activity within a 
clump, such as mass accretion, heats the surrounding dust caus- 



ing it to emit the 24-jt/m IR radiation. Towards eight clumps, 
a group of 8-/vm point sources is resolved. This indicates that 
the clumps were fragmented into smaller units, and the forma- 
tion of a stellar group/cluster is taking place. For comparison, 
from their total sample of 190 clumps associated with IRDCs, 
Chambers et al. (2009) found that 98 (52%) contain 24-jum 
sources, a percentage similar to ours. Parsons et al. (2009) found 
that from their (sub)sample of 69 clumps within IRDCs, 48 
(70%) were associated with embedded 24-jum source(s). More 
recently, Tackenberg et al. (2012) found that only ~ 23% of the 
clumps in their ATLASGAL-survey study showed no signs of 
IR emission. 

The estimated masses of the detected IR-dark and IR-bright 
clumps lie in the ranges ~ 25 - 2 x lO'* M© and ~ 42 - 9.5 x lO"' 
Mq, respectively. The clump masses as a function of radius are 
plotted in Fig.[TT] As can be seen, many of the clumps appear to 
lie above the mass-radius threshold for massive-star formation 
proposed by Kauffmann & Pillai (2010). In particular, 31 IR- 
dark clumps (~ 78%) lie on or above this threshold line, hence 
being potential sites of future high-mass star formation. For ex- 
ample, Beuther et al. (2011) estimated that the initial mass of a 
clump has to be ~ 850 Mq if it is to form at least one 20 Mq 
star; similarly, a ~ 1 900 Mq clump is needed to form a 40 Mq 
star. In this regard, our sample contains 28 (22) IR-dark clumps, 
which could be able to give birth to a 20 (40) Mq star. 

Following Casoli et al. (1986; their Sect. 4.1), we can com- 
pute the IR luminosities, Lir, of the detected IRAS sources using 
the sources' distances and flux densities at 12, 25, 60, and 100 
fim. The obtained Ljr values of IRAS 17474-2704, 18 112- 1720, 
181 14-1718, and 181 17-1738 are 4.4 x lO'*, 1.3 x lO'*, 7.4 x 10^ 
and 1.9 x 10^ Lq, respectively. These should be interpreted as 
lower limits to the bolometric luminosity, Lbob because other 
wavelength data (e.g., near-IR and submm) were not employed. 
These high luminosities are indicative of massive-star formation 
in the corresponding clumps. The clump G2. 1 1-SMM 5, associ- 
ated with IRAS 1747 4-2704, also shows C lass ii CH3OH maser 
emission at 6.7 GHz dCaswell et al. 19931 . Because this maser 
transition results from radiatively pumped population inversion 
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Fig. 9. Distributions of a) clump effective radii, b) masses, c) H2 column densities, and d) H2 number densities. Open histograms 
represent the IR-dark clumps, whereas shaded histograms represent the IR-bright clumps (those with 8/24 fim emission). The solid 
and dashed vertical lines indicate the average values for IR-dark and IR-bright clumps, respectively. The average (median) radii 
for IR-dark and IR-bright clumps are 1.19 (1.14) and 0.96 (0.87) pc, respectively. The average (median) mass of IR-dark clumps 
is 3 079 (2 108) Mq, whereas that of IR-bright clumps is 1 604 (804) Mq. The mean H2 column density is 3.2 x 10^^ cm"^ for 
both distributions (median values are 2.8 x 10^^ and 2.1 x 10^^ cm"^ for IR-dark and IR-bright clumps, respectively). The average 
(median) H2 number densities for IR-dark and IR-bright clumps are 7 250 (6 000) cm~^ and 6 667 (6 000) cm~^, respectively. 



by IR emission from warm dust, it can only occur near MYSOs 



(e.g., |Cragg et al. 1992[lMrnier et al. 20031 IXu et al. 2008^ . 

Krumholz & McKee (2008) showed that a clump's mass sur- 
face density needs to be E > 1 g cm"^ if heating is to prevent 
clump fragmentation into smaller units, and hence, enable high- 
mass star formation. In terms of H2 column density, this thresh- 
old is A^(H2) = E/jUh,»Jh - 2.1 x 10^^ cm"^. The estimated 
beam-averaged A^(H2) values are all smaller than this, but the 
(possible) substructure within the clumps, or dense cores, can 
have much higher column densities. Also, as discussed above, 
some of the clumps already show clear signposts of high-mass 
star formation. 



jChambers et al. 20091 1. The enhanced/extended 4.5-fj.m emis- 
sion is believed to be mainly caused by shock-excited H2 lines 
and/or ro-vibrational lines of CO(v - 1-0), implying the 
presence of outflow shocks in the source (IMarston e t al. 2004] 
Smith & Rosen 2005 ; Smith et al. 20061 lYbarra & Lada 200^91 
De Buizer & Vacca 2010). However, Takami et al. (2012) dis- 
cussed the possibility that the dominant emission mechanism 
responsible for EGOs might be scattered continuum in outflow 
cavities. In any case, EGOs appear to be related to outflow activ- 
ity. Finding EGOs among our sources is unsurprising, because 
they often appear to be associated with IRDCs and 6.7-GHz 
Class 11 CH3OH masers Chen etal. 20I0I I. 



5.1.1. Are there any EGOs ? 



5.1.2. The bubbles N10 and Nil 



A visual inspection of the Spitzer 4.5-fj.m images revealed 
an extended-like emission associated with the clumps SMM 
13, 17, 23, 25, 27, 29, and 32 in G13.22 (see Fig. [13. 
These type of sources are known as Extended Green 
Objects (EGOs; Cyganowski et al. 2008 1 or "green fuzzies" 



A zoom-in view towards the N 10/11 bubble environment in 
G13.22 is shown in Fig. [13] Our LABOCA clumps SMM 4, 5, 
7, and 1 1 appear to form a 7'4 or ~ 9.1 pc long "ridge" between 
these two IR bubbles from Churchwell et al. (2006). Moreover, 
the clumps SMM 6 and 9 form a filamentary structure, extend- 
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Right ascension (J2D0O) Right ascension (J2000) 




Riglit ascension (J2D00) 

Fig. 12. Spitzer-IRAC three-colour composite images towards EGO candidates in G13.22 overlaid with contours of LABOCA dust 
continuum emission. The 3.6, 4.5, and 8.0 jim emission is coded in blue, green, and red, respectively, and the colours are shown in 
linear scales. The contours are as in Fig. 2] Towards G13.22-SMM 13 {top left), the black circles indicate the positions of the 1.5- 
GHz VLA radio sources from Garwood et al. (1988), the diamond symbol shows the position of the 5-GHz VLA radio source from 
Becker et al. (1994), and the white cross shows the nominal catalogue position of IRAS 18114-1718. Towards G13.22-SMM 29 
(middle right), the position of IRAS 181 17-1738 is indicated by a black cross. The black circle towards G13.22-SMM 32 (bottom) 
denotes the 1.4-GHz radio source 181449-173243 from Condon et al. (1998). The white plus signs in the middle left and bottom 
panels represent the positions of our C'^O observations. We note that in the bottom panel the linear scale bar corresponds to the 
line-of-sight distance of SMM 25. 
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Fig. 11. Relation between mass and effective radius for the de- 
tected clumps. The IR-dark and IR-bright clumps have been 
plotted with open and filled circles, respectively. The dotted 
line represents the mass-radius threshold for massive-star forma- 
tion proposed by Kauffmann & Pillai (2010), i.e., MiR) = 870 



ing perpendicularly to the southeast from the above mentioned 
ridge. 

Churchwell et al. (2006) classified NIO and Nil as com- 
plete (or closed) bubbles, with NIO enclosing a star clus- 
ter (see also Watson et al. 2008, their Fig. 14). Moreover, the 
N 10/11 system was classified as a bipolar bubble whose lobes 
are in contact. Both bubbles are coincident with Hii regions 
(Deharveng et al. 2010| and references therein), and Watson et 
al. (2008) identified four possible ionising stars located inside 
NIO in projection. Watson et al. (2008) also identified four em- 
bedded candidate MYSOs on the rims of NIO, and suggested 
that the bubble could be associated with triggered massive-star 
formation. The elongated dust emission morphology observed 
here is perhaps consistent with this scenario. The third observa- 
tion noted by Watson et al. (2008) was that inside the 8-jum shell 
of NIO, emission at both 24 fim and 20 cm peak at the same 
position, implying the presence of hot dust inside the Hii region. 

Deharveng et al. (2010) reported the detection of two 
LABOCA condensations on the border of NIO. As shown in 
Fig. [13] the two condensations are coincident with our clumps 
SMM 5 and 7. As discussed by these authors, a Class ii methanol 
maser associated with the edge of SMM 5 (their condensation 
1) supports the scenario of triggered massive-star formation. 
Another CH3OH maser is seen close to the centre of NIO in pro- 
jection ( Pandian et al. 200^ . Deharveng et al. (2010) also spec- 
ulated that NIO could be in the process of opening, because of 
its elongated shape. They speculated that the observed bipolar 
morphology of the system might result from the expansion of an 
Hii region simultaneously in two opposite directions through the 
edge of the cloud. The dust filament consisting of SMM 6 and 
9 could perhaps be related to the bubble expansion. As already 
noted by Deharveng et al. (2010), Nil is not associated with 
significant LABOCA 870-/im emission. 



5.2. Lifetime of massive IR-dark clumps 

We can use the relative numbers of IR-dark and IR-bright clumps 
to estimate the statistical lifetime of the former stage. We avoid 
to use the term "starless" clump here, because some of the IR- 
dark clumps may well harbour faint YSOs, which cannot be de- 
tected with the current detection limit. 



Following Chambers et al. (2009), we also adopt as a re- 
presentative YSO lifetime the accretion timescale ~ 2 x 10^ yr 
from Zinnecker & Yorke (2007). If we further assume that /) 
all our IR-bright clumps, including those which are only partly 
associated with diffuse-like emission, host YSOs, ii) the star- 
formation rate is constant as a function of time, and ///) the 
clump's lifetime does not depend on its mass, we estimate the 
duration of the IR-dark phase of clump evolution to be TiR_daik ~ 
40/51 X 2 X 10^ ~ 1.6 X 10^ yr Conserning the third assump- 
tion here, Clark et al. (2007) noted that the clumps of different 
mass are expected to have different lifetimes, because their free- 
fall timescales, Tff oc 1 / ^Jp, can be different. They pointed out 
that this timescale problem may not be a problem if the studied 
clumps have roughly the same density, i.e., the same Tff. The es- 
timated average densities of our IR-dark and IR-bright clumps 
vary by factors of ~ 7 and 15, respectively, so their evolution- 
ary rates can also vary. This makes our third assumption rather a 
rough one. 

For comparison. Chambers et al. (2009), by considering only 
their quiescent (no embedded IR emission) and active (green 
fuzzies and 24 fim emission) sources (69 and 37 sources, re- 
spectively), estimated a factor of ~ 2 longer duration for the 
quiescent phase, i.e., ~ 3.7 x 10^ yr However, by comparing the 
relative numbers of their quiescent and 24-fim sources only ('in- 
termediate cores'), 69/98 - 0.7, the quiescent-phase timescale 
becomes ~ 1.4 x 10'' yr, which is very similar to our estimate. 
Parsons et al. (2009), adopting the embedded YSO duration of 
lO'* - 10^ yr, deduced that the IR-dark phase of massive clumps 
lasts a few times 10^ - 10"* yr With the aid of SiO spectral-line 
data, Russeil et al. (2010) estimated that the combined lifetime 
of starless clumps (no IR emission/high-velocity SiO emission) 
and IR-quiet MYSOs is ~ 1 x 10"* -H 6 x lO'* = 7 x lO'* yr in the 
NGC 6334/6357 complex, which is a factor of ~ 2 less than our 
T"iR-dark valuc. Tackcnbcrg et al. (2012) inferred the lifetime of 
massive IR-dark ('starless' in their nomenclature) clumps to be 
~ (6 ± 5) X 10"* yr, the upper limit being roughly comparable to 
our estimate. We also note that in their He rschel/Hi-GAL study 
of IRDCs, Wilcock et al. (2012) inferred the lifetime ~ 2 x 10^ 
yr for the IR-quiet phase, which is close to our estimate. We note 
that Wilcock et al. (2012) adopted the same YSO lifetime as we 
( ~ 2 X 10^ yr). 

Inspecting the diverse and somewhat heterogeneous lifetime 
estimates for IR-dark clumps, it seems that the true value is 
somewhere between ~ 10"^ - 10^ yr, although values as low as 
< 10^ yr have also been proposed (IMotte et al. 2007) . Extremely 
short lifetime of massive IR-dark clumps would, however, mean 
that it would be rather unlikely to see them in a large number in 
observational surveys. 



5.3. Clump mass distribution 

The clump mass distributions (CMDs) for the entire sample 
(composed of both IR-dark and IR-bright clumps), and sepa- 
rately for the IR-dark and IR-bright clumps are shown in Fig. [T4l 
The CMDs are plotted as dA^/dlogM versus M, where the 
first term is approximated as the number of clumps in each 
bin divided by the logarithmic mass interval, i.e., AN/AlogM. 
Following Lopez et al. (2011), we kept the histogram bin size 
A log M at a constant value of about 0.44. The CMDs were fitted 
with power-laws of the form dA^/d log M oc M^^, with the slopes 
r = 0.8 + 0.1, 0.7 ± 0.2, and 0.7 + 0.1 for the entire sample, IR- 
dark clumps, and IR-bright clumps, respectively. The errors of 
the slopes were determined by considering only the statistical 
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Fig. 13. Spitzer-IRAC three-colour composite image towards the bubbles NIO and Nil overlaid with contours of LABOCA dust 
continuum emission. The 3.6, 4.5, and 8.0 jim emission is coded in blue, green, and red, respectively, and the colours are shown in 
linear scales. The contours are as in Fig. |4] The squares show the positions of the two 870-yum condensations found by Deharveng 
et al. (2010); their condensation 1 is associated with our clump SMM 5, and condensation 2 is associated with SMM 7. The circles 
show the positions of the candidate embedded MYSOs on the rims of NIO (Watson et al. 2008 ). The diamond symbols indicate the 
positions of the 6.7-GHz Class ii methanol masers from Szymczak et al. (2002; right) and Pandian et al. (2008; left). 



Poisson uncertainty of the data, and the mass uncertain- 
ties were not taken into account. The mass ranges used in the 
fits were ~ 2.4 x 10^ - 1.8 x 10"* Mq for both the entire sam- 
ple and IR-dark clumps, and ~ 1.5 x lO"* - 1.1 x 10"* M© for 
the IR-bright clumps. Even the CMD slopes for the IR-dark and 
IR-bright clumps are similar, the CMDs themselves appear to 
be very different from each other; a two-sample Kolmogorov- 
Smirnov (K-S) test gives a probability of only 0.48% that they 
are drawn from the same underlying parent distribution. We note 
that the CMDs can also be expressed in the differential form of 
dN/dM oc M-", where a^F+l. 

The clumps studied here lie at different distances, and the 
sensitivities of the four maps are also different (median being 
~ 70 mJy beam"'). Therefore, the mass detection limit varies, 
and is not straightforward to determine. For example, at the me- 
dian distance of our sources, ~ 4.4 kpc, a source of 1 Jy cor- 
responds to a 15-K mass of ~ 225 Mq. Our completeness limit 
is, however, likely to be at a much higher mass, as suggested 
by the apparent peak near ~ 2 000 Mq. Because our sample is 
very heterogeneous in nature, it is not ideal for the CMD study. 
However, it is interesting to compare the derived slopes with 
those from the literature. 

The slope of the Salpeter (1955) initial mass function (IMF) 
is F = 1.35 or Q- = 2.35 for stars with masses in the range 0.4 
Mq < M < 10 Mq. Our slopes are clearly shallower than the 
Salpeter value. On the other hand, our results are similar to the 
CMDs of diffuse CO clumps, which are found to be well de- 
scribed by power-law forms with a between 1.6 and 1.8 (e.g., 
IStutzki & Glisten 19901 iKramer et al. 19981 ISimon et al. 2001l l. 



These slopes also resemble those of stellar clusters' mass func- 
tion, implying that massive clumps are not the direct progenitors 
of individual stars, but will instead presumably fragment to form 
groups/clusters of stars (e.g., |Elmegreen et al. 2000 



Rathbome et al. (2006) found a slope of or = 2.1 ± 0.4 above 
~ 100 Mq for their sample of MSX ^-fim dark clumps, which 
is quite close to the Salpeter power-law IMF, but also consistent 
with our values within the error bars. Interestingly, a two- sample 
K-S test yields practically a zero probability (~ 10"'~) that our 
IR-dark clump masses and those from Rathbome et al. (2006) 
are drawn from the common underlying distribution. We note 
that for this test the clump masses from Rathborne et al. (2006) 
were multiplied by 1 .2 14 to be consistent with the dust model we 
have adopted (their 1.2-mm dust opacity of 1.0 cm^ g"' was re- 
placed by the value /fi.2mm - 0.82 cm^ g ')■ Ragan et al. (2009) 
also built a CMD for their entire sample of IR-dark and IR- 
bright clumps (cf. left panel of our Fig. [14), and derived the slope 
Of = 1.76 ± 0.05 (from 30 to 3 000 Mq), which is very similar 
to our corresponding value of 1 .8 ± 0. 1 . Tackenberg et al. (2012) 
found a Salpeter-like slope of o- = 2.2 for the CMD of candidate 
massive starless clumps. We note that Salpeter-like slopes have 
also been determined for the mass distributions of more evolved 
massive clumps (e.g., IReid & Wilson 20051 IBeltran et al. 20061 
Bally et al. 2010 1, although the uncertainties are often too large 
to say whether they are actually flatter (or steeper) than the ex- 
act Salpeter slope. Indeed, in addition to our study, some au- 
thors have found the mass functions of evolved massive clumps 
to be similar to those of CO clumps (e.g., IBeuther et al. 20111 
Lopez et al. 2011 1. 
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Mass distributions of whole IRDCs have also been studied 
by Simon et al. (2006), Marshall et al. (2009), and Peretto & 
Fuller (2010), who determined the slopes a = 1.97 ± 0.09, 
a = 1.75 + 0.06, and a - 1.85 + 0.07, respectively. These are 
also comparable to our results and to the CMD slopes of the 
CO clumps. Marshall et al. (2009) suggested that the similar- 
ity between the mass distributions of IRDCs and CO clumps 
perhaps indicates that IRDCs are the result of density fluctua- 
tions caused by interstellar turbulence. The analytical theory by 
Hennebelle & Chabrier (2008) is able to explain the flat mass 
spectra of CO clumps formed by supersonic turbulent flows. In 
the case of large-scale turbulent flows, the slope of the CMD's 
power-law tail is determined by the spectral index of the turbu- 
lent power spectrum. In our notation, the CMD slope F in the 
Hennebelle-Chabrier theory is given by F = 2 - n'/3, where «• is 
the three-dimensional power-spectrum index of the logarithmic 
density field. As discussed by Hennebelle & Chabrier (2008), 
«• appears to be ^ 11/3, i.e., similar to the Kolmogorov power 
spectrum index 11/3 of the velocity field in incompressible tur- 
bulence. This leads to the value F ^ 7/9 ^ 0.8, which is re- 
markably close to the CMD slopes of CO clumps and the values 
derived in the present study. 



5.4. On the origin and iiierarciiical fragmentation of 
filamentary IRDCs 

The formation of filamentary structures by colliding shock 
fronts is clearly seen in numerical simulations (e.g.. 



Banerjee et al. 2006i. Large-scale colliding supersonic tur- 
bulent flows offer an intriguing mechanism for the formation of 
filamentary IRDCs. As discussed in the previous section, the 
IRDCs' mass functions can be understood in terms of turbulent 
flows. Moreover, Jimenez-Serra et al. (2010) found extended 
SiO emission along the filamentary IRDC G035. 39-00.33, 
which might have been produced in shocks caused by converg- 
ing flows. Hernandez & Tan (2011) and M12 found evidence 
that the surface pressure may play an important role in the 
dynamics of filamentary IRDCs. This could also be related to 
the cloud formation in gas overdensities resulting from colUding 
flows. 

Although the formation of IRDCs could be caused by in- 
terstellar turbulent flows, their further fragmentation may have 
its origin in some other mechanism(s). As was shown in the 
present study, the fragmentation of the filamentary IRDC Gl 1 .36 
into clumps can be explained by the "sausage"-type fluid in- 
stability. In an IRDC, this was first found to be the case in 
the Nessie Nebula by Jackson et al. (2010). More recently, 
M12 inferred that this is also the case in the filamentary IRDC 
G304. 74-1-01. 32. Hence, to our knowledge there are three dif- 
ferent filamentary IRDCs where the predictions of the sausage- 
instability theory have been tested so far. In all cases, the theory 
and observations are in very good agreement with each other 
One remarkable feature that should be raised is that in all the 
three above mentioned cases the fragmentaion length scale cor- 
responds to the wavelength of the fastest growing mode only 
when the non-thermal motions are taken into account. Non- 
thermal ('turbulent') motions are therefore still important in the 
fragmentation process, on least at the scale of clumps. However, 
this can also continue in a hierarchical way down to the scale of 
cores, as found by Wang et al. (201 1) in the massive IRDC clump 
G28. 344-0. 06-Pl . A Jeans-type gravitational instability may start 
to dominate inside dense cores, i.e., cause them to fragment into 
still smaller condensations (cf . IMiettinen et al. 2012l l. 



5.5. CO depletion and its implications on the age of the 
G11. 36+0.80 filament 

The CO depletion factors we estimate towards 25 target posi- 
tions lie in the range ~ 1.0 + 0.2 - 20.7 ± 0.8. Some of the val- 
ues are very large, taking that the beam size 27'.'8 is probing the 
scale of clumps at the source distances. In contrast, using the 
same kind of APEX/C'^O data as here, Miettinen et al. (2011) 
found depletion factors of only ~ 0.6 - 2.7 towards a sample of 
seven clumps in IRDCs and similar values of ~ 0.3 - 2.3 were 
found by M12 towards the clumps in the IRDC G304.74. On the 
other hand, Chen et al. (201 1), who used 34" C'^O observations, 
reported /d values as high as ~ 19 in the IRDC G34.43-i-0.24, 
which resembles our largest /d values. Pillai et al. (2007) de- 
rived the CO depletion factor of 9.3 towards one of our clumps, 
namely G13.22-SMM 5 (their source G13. 18+00.06). 

Some of the largest /d values in the present study are 
found towards the positions C, D, E, and F in the filamentary 
IRDC in G1.87, and position A in G13.22. These might be 
potential sites to search for high levels of molecular deutera- 
tion, since the depletion of CO plays a key role in the deu- 
terium fractionation. Interestingly, the largest fu value of ~ 
20.7 is seen towards the edge of G2.11-SMM 5 (position H). 
This raises the question whether the high depletion factor is 
the result of molecular freeze-out onto dust grains, or pos- 
sibly due to photodissociation by FUV photons at the edge 
of the clump (ITielens & Hollenbach 19851 IVisser et al. 2009] l ? 
Photodissociation could also play a role in reducing the C'^O 
abundance in some other positions of our line observations. 

The fu values in the Gl 1.36 filament (excluding positions F 
and G), where we best probe the submm peak positions among 
our sources, are ~ 1.6 + 0.2 - 3.9 ± 0.6. This implies that non- 
negligible CO depletion might be present in some parts of the 
filament. Hernandez et al. (2011) inferred depletion factors of 
~ 3 - 4 for the thinnest part of the filamentary IRDC G035.30- 
00.33, and used the CO depletion timescale to estimate the 
cloud's age. Here, we do the same exercise for Gl 1.36. The CO 
depletion timescale can be computed as the inverse of the freez- 
ing rate ( [RawUngs et al. 1992| as 



1 



1 



^dep — 



^freeze ngTTa^VcO^ 



(3) 



where Kg is the grain number density, Og is the average grain 
radius (na^ is the mean cross section of the dust grains, and 

rtgTTal is the grain opacity), vco is the thermal speed of the CO 
molecules, and S is their sticking coefficient. We can write the 
grain number density as rig = Xgn(H2), where Xg is the fractional 
grain abundance given by 



WH2 „ '"Hi^d 

Xg = xRi - 

m„ 



3 "g/^g 



(4) 



In the above formula, niuj is the H2 molecule's mass, nig (pg) 
is the mass (density) of the grain particle, and = 1/100 
is again the dust-to-gas mass ratio. By adopting the "stan- 
dard" values flg = 0.1 fim and pg - 3.0 g cm"-', we obtain 
Xg - 2.7 X 10"'^. The Maxwellian mean thermal speed of 
the CO molecules (mco = 28mH) at Tki,, - 15 K is vco = 
(8feB7'kin/7rmco)'^^ - 106 m s"'. Assuming that 5 = 1, i.e., 
that the CO molecules stick to the dust grains in each collision, 
we can write Tjep ~ 4.4 x 10''/n(H2) [cm"-'] vf^ The volume- 



We note that the formula used by Hernandez et al. (2011; their 
Sect. 4.2) gives a factor of 5.5 times shorter CO depletion timescale. 
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Fig. 14. From left to right are shown the differential (dA^/dlogM) clump mass distributions of the entire sample, IR-dark clumps, 
and IR-bright clumps, respectively. In all panels the mass bin size is A log(M/MQ) ^ 0.44. The error bars correspond to the standard 
Poisson Va^ counting uncertainties. The solid lines indicate the best- fit power-law functions of the form dN/ d log M oc M^^, where 
the slope is given in each panel. The fits refer to masses above ~ 2400 Mq for both the entire sample and IR-dark clumps, and 
above ~ 1 500 Mq for the IR-bright clumps. We note that the Salpeter IMF has a slope of F = 1.35. 



average density in the Gil. 36 clumps is ~ 10"^ cm ^ (Table [3]), 
which gives Tdep ~ 4.4 x 10^ yr Interestingly, this is very close to 
the estimated fragmentation timescale in Gil. 36 (~ 5.4x10^ yr). 
The age of the filament may therefore be similar to its fragmen- 
tation time. In this case, the cloud's fragmentation would have 
started at the same time as it was formed. Perhaps the cloud- 
forming flows also triggered the sausage-type fluid instability in 
the cloud. 



6. Summary and conclusions 

Four selected regions in the Galactic plane, all contain- 
ing IRDCs, were mapped with APEX/LABOCA at 870 fim. 
Moreover, selected positions in the fields were observed in 
C'^0(2 - 1) with APEX. Our main results and conclusions can 
be summarised as follows: 

1. The total number of clumps identified in this survey is 91. 
The number of Spitzer 8- and lA-jiva dark clumps was esti- 
mated to be 40. The remaining 5 1 clumps were found to be 
associated with both 8 and 24 jim emission, or only with 
24 fiva emission (embedded point sources, group of point 
sources, diffuse/extended emission). 

2. Many of the clumps appear to be potential sites of (future) 
high-mass star formation, as implied by the comparison with 
the mass-radius threshold proposed by Kauffmann & Pillai 
(2010). Some of the clumps already show clear signposts of 
ongoing high-mass star formation, such as Class ii CH3OH 
maser emission. 

3. Seven clumps show extended-like 4.5 /im emission, hence 
are classified as candidate EGOs. 

4. One of the mapped fields (G13.22) contains the Spitzer MIR- 
bubbles NlO/11 found by ChurchweU et al. (2006). The 
bubble region is believed to be associated with triggered 
massive-star formation, and the observed dust emission mor- 
phology might have been created by the bubble shell-cloud 
interaction process. 

5. The relative numbers of IR-dark and IR-bright clumps im- 
plies the duration of the IR-dark stage to be ~ 1.6 x 10^ yr 
Although this estimate is likely to suffer from the relatively 
low number statistics, it is comparable with some earlier life- 
time estimates ^Chambers et al. 20091 IWilcock et al. 2012] |. 
Surveys with higher number of sources typically indi- 
cate shorter Ufetimes for massive IR-dark clumps (e.g., 
Tackenberg et al. 2012 1. 



6. The clump mass distributions were constructed for the to- 
tal clump sample, and for the IR-dark and IR-bright clumps 
separately. In all cases, the high-mass tail could be fitted with 
the power-law function of the form dA/d log M oc M^^ , with 
r ^ 0.7 . . .0.8. This is very similar to the mass functions 
found for the diffuse CO clumps, and can be understood 
in terms of supersonic-turbulence induced cloud formation 
jHennebeUe & Chabrier 2008l l. 

7. The C'^O observations revealed potential targets of strong 
CO depletion. In some cases, however, our line observations 
probe more dilute gas, where the low C'^O abundances could 
be caused by photodissociation by FUV photons. 

8. The filamentary IRDCG11.36H-0. 80 was studied in more de- 
tail, because there our line observation positions were well 
matched with the submm peak positions. The filament's frag- 
mentation into clumps can be well explained by a "sausage"- 
type fluid instability. In particular, the observed projected 
clump separations are in excellent agreement with the the- 
oretical prediction. This is in agreement with the results for 
other IRDCs (Jacks on etal. 20101 M12). Most clumps in the 
filament appear to be gravitationally bound, and the fila- 
ment as a whole appers to be close to virial equilibrium. 
Interestingly, the estimated fragmentation timescale of the 
filament and the CO depletion timescale were inferred to be 
comparable to each other (~ 5 x 10^' yr). The perturbations 
responsible for the filament fragmentation into clumps might 
have been excited by the cloud-forming process, i.e., by con- 
verging turbulent flows. 
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Table 3. Characteristics of the clumps identified from the LABOCA maps. 
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point(fg)^ /dark 


SMM 40 . . . 


17 50 48.7 


-27 22 01 


0.52 ±0.10 


0.43 ± 0.07 


12.2 


0.63 ± 0.03 


356 ± 67 


1,6 ± 0,3 


0,7 ±0,1 


SSTGLMC GOOl.9357-00.1696' (4'.'1); 
SSTGLMC GOOl.9323-00.1700' (8'.'7); 
BGPS Gl.934-0.1722 (ia'2) 


group/point 


G2.11+0.00 
























SMM 1 . . . 


17 50 26.8 


-27 02 16 


0.37 ± 0.07 


2.41 ±0.25 


36.8 


1.11 ±0.09 


846 ±291 


1.8 ±0.3 


0.3 ±0.1 




dark/dark 


SMM 2(D)'' . . . 


17 50 30.6 


-27 07 20 


0.30 ± 0.07 


0.28 ± 0.07 


12.4 


0.37 ± 0.03 


105 ± 42 


1.5 ±0.3 


1.0 ±0.4 




dark/dark 


SMM 3 . . . 


17 50 30.9 


-27 04 12 


0.47 ± 0.08 


3.17 ±0.32 


37.2 


1.12 ±0.09 


1 113±381 


2.3 ± 0.4 


0.4 ±0.1 


BGPS G2. 154+0.036- (6'.' 2) 


dark/dark 


SMM 4... 


17 50 32.1 


-27 08 48 


0.37 ± 0.07 


1.50 ±0.16 


29.0 


0.87 ± 0.07 


527 ± 181 


1.8 ±0.3 


0.4 ±0.1 


BGPS G2.090-0.0082 (16'.'0) 


point{fg)^ /dark 


SMM 5... 


17 50 36.0 


-27 05 44 


2.74 ± 0.28 


6.27 ± 0.63 


34.4 


1.23 ± 0.02 


1 457 ± 152 


4.9 ± 0.5 


0,4 ±0,1 


GPSR5 2.1434-0.010" (1'.'6); 
G002.143+00.0105(2'.'5); 


extended/point 



G002.14+0.01 A*(3'.'0); 
IRAS 17474-2704^ (3'.'7); 



Table 3. continued. 



Source 


azooo.o 


^2000.0 


'S70;»m 




Relf 


M 


N(H2) 


(n(H2)> 


Associations" 


8/24 ijm 




[h:m:s] 


[":':"] 


[Jybeam"'] 


[Jy] 


["] [pc] 


[Mq] 


[10^^ 


[10" cm-3] 







2MASS J17503612-2705483 (4'.' 4); 

002.14+00.01* (5'.'9); 
MSX5C 6002.1419+00.0099' (6'.' 6); 
BGPS G2.144+0.0062 (13'.'8); 
002.14+00.01"' (16'.'0) 



SMM 6 ... 


17 50 36.9 


-27 06 32 


0.38 ± 0.07 


1.69 ±0.18 


28.8 


0.87 ± 0.07 


379 ± 131 


1.2 ± 0.2 


0.3 ±0.1 


■■■ 


point/point 


SMM 7 ... 


17 50 41.7 


-26 59 28 


0.33 ± 0.07 


1.46 ±0.16 


28.8 


0.87 ± 0.07 


513 ± 177 


1.6 ± 0.3 


0.4 ±0.1 


BGPS G2.240+0.046^ (6'.'2) 


point(fg)^ /dark 


SMM 8 . . . 


17 5103.3 


-27 11 24 


0.39 ± 0.07 


1.47 ±0.16 


26.8 


0.81 ± 0.06 


516 ± 178 


1.9 ± 0.3 


0.4 ± 0.2 


BGPS G2.1 14-0.130^ (11'.'2) 


dark/dark 


SMM 9 ... 


17 51 15.3 


-27 03 19 


0.32 ± 0.07 


1.64 ±0.17 


30.4 


0.91 ± 0.07 


576 ± 198 


1.6 ± 0.3 


0.4 ±0.1 




weak/diffuse 


SMM 10 . . . 


17 51 16.5 


-27 05 39 


0.38 ± 0.07 


2.57 ± 0.26 


36.6 


1.10 ± 0.09 


902 ± 309 


1.9 ± 0.3 


0.3 ± 0.1 




dark/dark 


(jll.36+U.»U 
























oMM 1 (A) . . . 


18 U/ 5d.\j 


-18 43 46 


U.6z ± U.U/ 


1.58 + 0.16 


23.9 


A 1 A A/T 

0.3s ± O.Oo 


128 ± 41 


2.0 ± 0.2 


1.1 ± 0.3 




point/point 


SMM 2{C f . . . 


lo 07 35.6 


10 11 

-Is 43 22 


u.4s ± u.u6 


11/1 1 A 1 1 

1.24 ±0.13 


23.0 


0.37 ± 0.06 


157 ± 50 


2.4 ± 0.3 


1.4 ± 0.5 




dark/dark 


SMM J(B) . . . 


18 07 35.8 


10 A'^ A'^ 

- 1 8 42 42 


0.75 ± o.oy 


2.28 ± 0.23 


29.6 


0.47 ± 0.07 


184 + 59 


2.4 ± 0.3 


0.8 ± 0.3 


sue Cil 1.360+0.800 (3. 1) 


point/point 


SMM 4 . . . 


18 07 36.1 


to A A 

- 1 8 44 26 


0.60 ± 0.07 


1.25 ±0.13 


23.5 


0.38 ± 0.06 


101 ±32 


1.9 ± 0.2 


0.8 + 0.3 


SSTGLMC GOl 1 .3364+00.786P (2'.'2) 


group/group 


SMM 5(5/' . . . 


18 07 36.7 


-18 41 14 


0.62 + 0.07 


2.37 ± 0.24 


32.5 


0.52 ± 0.08 


192 ± 61 


2.0 ± 0.2 


0.6 ± 0.2 


SDC Gil. 380+0. 809'^ (lO'.'l) 


point/point 


SMM 6(¥f . . . 


18 07 39.0 


-18 42 10 


0.20 ± 0.04 


0.15 ± 0.04 


11.1 


0.18 ± 0.03 


25 ± 9 


1.0 ± 0.2 


2.0 ± 0.7 


SDC Gl 1.374+0.792^^ {9'.'1) 


dark/dark 


SMM 7(Gr . . . 


18 07 40.4 


-18 43 18 


0.36 ± 0.05 


A A AO 

0.72 ± 0.08 


21.2 


f\ A .A f\C 

0.34 ± 0.05 


58 ± 19 


1.1 ± 0.2 


0.7 ± 0.2 


SDC Gll.361+0.777 (12. 8) 


poinVpoint 


i7lO.2Z-U.UO 
























SMM 1 . . . 


1 1 "3 An r\ 
18 13 4/.U 


-1 / 11 09 


0.49 ± 0.09 


A jCI 1 A 1 f\ 

0.62 ± 0.10 


15.2 


A 1 A AT 

0.32 ± 0.02 


85 ± 18 


1.5 ± 0.3 


1.2 ± 0.3 


BCjFS G13.245+0.158 (5. Oj 


poinVpoint 


SMM 2 ... 


18 13 53.4 


-17 31 10 


f\ A A .A AA 

0.44 ±0.09 


f\ C A A 1 A 

0.54 ± 0.10 


14.8 


A '^A . A A1 

0.30 ± 0.02 


72 ± 17 


1.4 ± 0.3 


1.2 ± 0.3 


JCMTSF J181353. 2-1731 19 (14. 7); 


dilfuse/weak 






















SDL G13. 125+0. 072 (2(7. 5) 




oMM D . . . 


18 13 54.6 


1 "7 OA /to 


A /I /I _i_ n no 


1 AO _i_ A 1 /I 

1 .uy ± u. 14 


21.4 


A /I ^ _i_ A A'J 
U.4j ± U.U3 


150 + 28 


1 .4 + 0.3 


0.8 + 0. 1 




■ 

poinl/point 


>3 iVI IVI 4 . . . 


18 13 jj.y 


1 1 OQ lA 

-1 / 28 34 


1.21 ± 0. 15 


y.i2 ± u.y2 


45.2 


A Q'2 J. A an 
U.yj ± 0.0/ 


1 JlU ± ZZZ 


3.8 ± 0.5 


0.1 ±0.1 


J CM 1 or J 181356. 1 -1 1 Lai 1 (4. 3 ), 


poinl/point 






















MIA /I 1 2 /c" 1 \. 

IN lU-4 \j. I}, 
























MIA "J 12 (I" Q\ 
IN lW-3 \i -o} 




SMM 5 . . . 


1 1 /I c\c\ n 
18 14 OU. / 


-1 / z8 38 


5.61 ± 0.57 


25.92 ± 2.59 


52.7 


1 AO ■ A AO 

1.08 ± 0.08 


5 335 ± /04 


17.0 ± 1.7 


1.2 ± 0.3 


CST\i~^ 1 '! 1 11 1 A A/; ill /ttt 1 \. 

SDL G13.1 / /+0.061 (2. Ij; 


point/extended 






















outilow (4. /); 
























C/^ A * jfTlC 1 1 O 1 A A*; 14 fCllfW, 

SLAMFS G13. 18+0.06 (5. 0); 
























JLMlSr J 181400. /-1 72843 (5. D), 
























INIO-/ (5. 6j; 
























G013.177+00.059SMM^^ (9. 1); 
























BGPS G13. 179+0.060^ (lft'5); 
























SSTGLMC G013. 1818+00. 0610' /NlO-9'^ (17'.'4) 




SMM 6 . . . 


18 14 08.8 


-17 28 57 


A O'^ i A 1 ^ 

0.83 ± 0.12 


A OA A CA 

4.89 ± 0.50 


36.2 


r\ n A .A AC 

0.74 ± 0.05 


652 ±114 


2.5 ± 0.4 


0.7 ±0.1 


T> T>o lAi , r\ r\'^ aI / A// <y\. 

BGPS G13. 191+0.034 (4. 8); 


extended/extended 






















JCMTSF J181408. 7-172907' (9'.'6); 
























G013. 19+00.04 b^^ {W.'A) 




SMM 7 . . . 


18 14 09.4 


-17 27 21 


3.02 ± 0.31 


48.22 ± 4.82 


75.2 


1.55 ± 0.11 


6433 ± 1 114 


9.5 ± 1.0 


0.8 ±0.1 


2MASS J18140960-1727219'^ (3'.' 6); 


extended/extended 






















JCMTSF J181409. 5-172725' (5'.'0) 




SMM 8 . . . 


18 14 09.9 


-17 19 53 


0.53 ± 0.10 


0.41 ± 0.09 


11.9 


0.25 ± 0.02 


42 ± 10 


1.8 ± 0.3 


1.2 ±0.3 


IRAS 18112-1720 (5'.'1); 


extended/extended 






















GPSR5 13.322+0.096'^ (8''2) 




SMM 9 . . . 


18 14 10.5 


-17 29 50 


0.78 ±0.11 


2.12 ± 0.23 


26.1 


0.54 ± 0.04 


295 ± 52 


2.5 ± 0.3 


0.9 ± 0.2 


SDC 013.177+0.017'^ (3''0); 


extended/extended 






















JCMTSF J181410.4-172955' (15'.' 4) 




oMIVl lU . . . 


12 ^A no 

io i'^ iz.z 


1 1 1 A 
-L / Lj 1.4 


U.j4 ± U.iU 


1 Qo J. n on 




n ^'5 J. r\ r\A 


yyz ± /u 


7 J. n ^ 
z. / ± U.J 


l.z ± u.z 


JLlVllor JIoI4IZ.D-i / ZD ZD (.1^. z). 


dark/dark 






















BGPS G13.253+0.044' (2a'2); 




SMM 11 ... 


18 14 14.4 


-17 26 30 


0.49 ± 0.09 


4.08 ± 0.42 


39.3 


0.81 ± 0.06 


544 ± 95 


1.5 ± 0.3 


0.5 ±0.1 




diffuse/diffuse 


SMM 12 . . . 


18 14 17.7 


-17 28 58 


0.45 ± 0.09 


1.99 ±0.21 


29.0 


0.60 ± 0.04 


265 ± 47 


1.4 ± 0.3 


0.6 ±0.1 


Unknown object" (4" 4); 


dark/extended 






















JCMTSFJ181417.9-172908' (10".' 6) 




SMM 13 . . . 


18 14 24.5 


-17 17 30 


1.41 ±0.16 


4.01 ±0.41 


29.0 


0.27 ±0.10 


65 ±47 


2.0 ± 0.2 


1.5 ±1.1 


GPSR5 13.385+0.066"* (4'.'8); 


extended/extended 






















013.385+0.069''' (5"2); 
























013.386+0.069''' (7'.' 6); 
























BGPS G13.387+0.0662 (7"8); 
























013.386+0.065'" (7'.'9); 
























IRAS 18114-1718 (ir.'O) 




SMM 14 . . . 


18 14 25.6 


-17 22 46 


0.51 ±0.09 


1.50 ±0.17 


23.0 


0.50 ± 0.03 


221 ± 39 


1.6 ± 0.3 


0.8 ±0.1 


SDC G13.310+0.021" (2'.'2); 


dark/point 


SMM 15 . . . 


18 14 26.7 


-17 34 58 


0.91 ±0.12 


2.58 ± 0.27 


28.2 


0.49 ± 0.05 


243 ± 60 


2.9 ± 0.4 


0.9 ± 0.2 


SSTGLMC G013.1343-00.07943 (2'.'5); 


poinVpoint 


SMM 16 . . . 


18 14 27.8 


-17 22 38 


0.54 ±0.10 


2.57 ± 0.27 


29.3 


0.63 ± 0.04 


379 ± 65 


1.7 ± 0.3 


0.7 ± 0.1 




poinVpoint 


SMM 17 . . . 


18 14 28.3 


-17 36 14 


1.23 ±0.15 


5.16 ±0.52 


33.6 


0.58 ± 0.07 


485 ± 120 


3.9 ± 0.5 


1.1 ±0.3 


SSTGLMC GO 13. 11 82-00.0966-^ (3'.'7); 


poinVpoint 



I 

i 

O 

n 
> 



a" 

ens 



n 



BGPS G13.121-0.094^ (19'.'0); 



Table 3. continued. 



Source 


^^"2000.0 


^2000.0 


870/(m 


■^870/(111 






M 


A'(H2) 


<«(H2)) 


Associations" 


8/24 /jm 




[h:m:s] 


[":':"] 


[Jy beam"^] 


[Jy] 


["] 


[pc] 


[Mq] 


[10^2 cm-2] 


[10* 


























SDCG13. 121-0.091" (19'.'9) 




SMM 18(A)''... 


18 14 28.6 


-17 33 26 


1.00 ±0.13 


19.29 ± 1.93 


49.2 


1.07 ±0.07 


2 860 ± 479 


3.1 ±0.4 


1.1 ±0.2 


SDCG13. 158-0.073" (5'.'6) 


extended/extended 


SMM 19... 


18 14 28.9 


-17 21 30 


0.40 ± 0.09 


1.08 ±0.13 


21.5 


0.46 ± 0.03 


249 ± 45 


2.0 ± 0.4 


1.2 ±0.2 




dark/dark 


SMM 20 . . . 


18 1431.4 


-17 36 06 


0.48 ± 0.09 


3.01 ±0.31 


33.4 


0.58 ± 0.06 


443 ± 109 


2.4 ± 0.4 


1.0 ±0.3 




dark/dark 


SMM 21 ... 


18 14 34.5 


-17 34 34 


0.48 ± 0.09 


2.09 ± 0.22 


28.2 


0.61 ± 0.04 


310±53 


1.5 ±0.3 


0.6 ±0.1 




point/point 


SMM 22(C)* . . . 


18 14 35.4 


-17 30 42 


0.55 ±0.10 


1.87 ±0.20 


24.5 


0.42 ± 0.05 


275 ± 68 


2.7 ± 0.5 


1.7 ±0.4 


BGPS G13.213-0.076- (5'.'9) 


dark/dark 


SMM 23(E)'' . . . 


18 14 36.8 


-17 29 22 


0.67 ±0.10 


1.90 ±0.21 


25.1 


0.43 ± 0.05 


179 ± 45 


2.1 ±0.3 


1.0 ±0.3 




point/point 


SMM lA ... 


18 14 37.0 


-17 38 50 


1.61 ±0.18 


8.35 ±0.70 


43.6 


2.62 ± 0.08 


9496± 1005 


5.1 ±0.6 


0.2 ±0.1 


SSTGLMC G013.0970-00.14472° (7'.'7); 
BGPS G13.()97-0.1462 (18'.'3) 


point/point 


SMM 25 . . . 


18 14 39.0 


-17 33 02 


1.69 ±0.19 


12.70 ± 1.27 


53.1 


1.15 ±0.08 


I 883 ±315 


5.3 ±0.6 


0.6 ±0.1 


BGPSG13.186-0.1082 (9'.'2); 
SDCG13. 190-0.105" (19'.'8) 


extended/point 


SMM 26 . . . 


18 14 39.6 


-17 41 10 


0.64 ±0.10 


2.51 ±0.26 


27.8 


1.67 ±0.05 


4 465 ±545 


3.2 ± 0.5 


0.4 ±0.1 


SDCG13.067-0.172"(4'.'4) 


dark/dark 


SMM 27(F)'' . . . 


18 14 40.7 


-17 29 22 


3.16 ±0.33 


13.25 ± 1.33 


44.3 


0.76 ± 0.09 


I 246 ± 307 


10.0 ± 1.0 


1.3 ±0.3 


SDCG13.246-0.081" (18'.'1); 
BGPS G13.245-0.0842 (19'.'9) 


group/point 


SMM 28 . . . 


18 1441.2 


-17 23 18 


1.53 ±0.17 


8.21 ±0.82 


38.8 


0.84 ± 0.06 


I 212 ±203 


4.8 ± 0.5 


0.9 ±0.2 


BGPS G13.333-0.0382 (16'.' 6) 


group/group 


SMM 29 . . . 


18 14 42.1 


-17 37 06 


2.42 ± 0.25 


11.27± 1.13 


46.3 


2.77 ± 0.09 


7 269 ±868 


8.3 ± 0.9 


0.2 ±0.1 


BGPSG13.133-0.1502(2'.'9); 
IRAS 18117-1738^' (11'.'2) 


extended/point 


SMM 30 . . . 


18 14 42.1 


-17 2145 


0.43 ± 0.09 


0.66 ±0.10 


16.3 


0.35 ± 0.02 


97 ±20 


1.4 ±0.3 


1.0 ±0.2 


SDC G13.357-0.029" (6'.'3); 
BGPS G13.359-0.0302 (13'.' 1) 


point/point 


SMM 31... 


18 14 46.0 


-17 23 53 


0.58 ±0.10 


1.87 ±0.20 


25.8 


0.56 ± 0.04 


432 ± 74 


2.1 ±0.4 


1.1 ±0.2 




dark/dark 


SMM 32'' . . . 


IS 14 49.9 


-17 32 45 


4.28 ± 0.44 


31.75±3.18 


62.3 


1.33 ±0.09 


2 620 ±443 


7.7 ±0.8 


0.5 ±0.1 


G0I3.2I0-0.14422(1'.'4); 


extended/extended 



SMM 33 . 
SMM 34 . 



18 14 52.4 
18 14 53.6 



-17 41 14 
-17 40 06 



0.55 ±0.10 
0.84 ±0.12 



0.98 ±0.13 
1.75 ±0.19 



18.2 
22.7 



1.09 ±0.04 
1.36 ±0.04 



I 743 ± 257 
1990 ±251 



2.7 ± 0.5 
2.6 ± 0.4 



0.6 ±0.1 
0.4 ±0.1 



MSX6C G013.2097-00.14362^ (1'.'6); 
PMNJ1814-17322'*(1'.'8); 
181154.9-173348^^ (5'.'4); 
013.210-0. 1452" (6'.'7); 
BGPSG13.211-0.1422(11'.'0) 
BGPS G13.091-0.2182(12'.'l); 
BGPSG13.109-0.2162(8'.'7) 



dark/dark 
point/point 



CD 
3 



I 

o 
o 

> 

3 



Notes. Remarks on the associated sources. Offset ["] from tfie 870-yum peak position is indicated in parenthesis.**' The clump is associated with the marked C'''0(2 - l)-observation ta^ 
get position.*''' The 8-yum point source is associated with SSTGLMC GOOl. 8134-00.1136, which has the [3.6] - [4.5], [4.5] - [5.8], and [4.5] - [8.0] colours of 0.098, 0.245, and 0.031, reS, 
spectively, and is likely to be a foreground star.'"" The 8-yum point source corresponds to SSTGLMC GOOl.9301-00.1234 with the colours [3.6] - [4.5] = 0.081, [4.5] - [5.8] = 0.158, an(p3 
[4.5] - [8.0] = 0.167, and is likely a foreground star.*'"' The 8-yum point source corresponds to SSTGLMC GOO 1.97 16-00. 1375 with the colours [3.6] - [4.5] = -0.033, [4.5] - [5.8] = 0.294O 
and [4.5] - [8.0] = 0.141, and is likely a foreground star.<^> The 8-yum point source is associated with SSTGLMC G002.0927-00.0039, which has the [3.6] - [4.5], [4.5] - [5.8], and [4.5] - [8.0? 
colours of 0.802, 0.775, and 0.290, respectively, and is likely to be a foreground star.<s> The 8-yum point source corresponds to SSTGLMC G002.2440-l-00.0447. Its [3.6] - [4.5], [4.5] - [5.8], and 
[4.5] - [8.0] colours are, respectively, 0.523, 0.436, and 0.838, and it is likely a foreground star.**' The clump can be resolved by eye into two "subclumps", but they are dealt as a one source 
by clumpfin d.<" SCUBA clump dPi Francesc"oeta l. 2008l.<^> Bolocam clump l |Rosolowsky et al. 2010| l.'^' YSO candidate iRobitaille et al. 2008t .<'^' UC Hll region JBecker et al. 1994> .'^' OH 
maser ( |Argon et al. 2000)."'' UC Hll region fForster & Caswell 2000 1."' MacLeod et al. (1998) search for 6.7-GHz CH3OH m aser emission towards this sources, but did not detect it.<^' 6.7- 
GHz CH^OH maser dCaswell et al. 1 995 1.'*" IR source (Egan et al. 2001 1.*"" OH and H2O maser source (see lCaswell et al. 1983>."" IRDC from the Peretto & Fuller (2009) catalogue; available 
at www.irdarkclouds.org/.*'^' YSO candidate dWatson et al. 2008i .*'^' CO outflow associated with the bubble NIO (Beau mont & Williams 2010i l.'''*' A submm clump studied by Pillai et al. 
(2007).*'^' SCUBA clump ([Thompson et al. 2006|l.<'*' Possible AG B star dCutnet al. 2003).*'^' Source nro. 28 in Russeil (2003); its nature is unknown.*"*' Radio source ( Becke r et al. 1994 !.<"' 1.5- 
GHz radio source (Garwood et al. 1988i.''^"* Possible AGB star dRobitaille et al. 20081 1.*^'' The distance ambiguity of the IRAS source was resolved by Sewilo et al. (2004).*-^' Hll region 
dWhite et al. 20051 [Urquhart et al. 2009; .'^^^' YSO ( ,Egan et al. 2003) .'^'" 1.4-GHz radio source dCondon et al. 1998< .'^^' Hll region tChini et al. 1987b .'^''' Hll region OWinketal. 1982t . 



